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Preface 


In writing this contribution to the Problem Solvers series, I 
have assumed that the reader has completed a first course in 
mechanics and is familiar with the techniques employed to 
solve problems relating to the motion of a single particle or a 
small number of interacting particles. In particular. I have 
supposed that his knowledge of the subject embraces inertial 
frames, Newton’s laws, the velocity and acceleration vectors 
and their components, properties of idealized components of 
mechanical systems such as inextensible strings, perfectly 
rough planes, etc., and the definitions of work and energy. 
However, some of the problems included in Chapter 2 serve 
to recall these fundamental concepts to mind. 

The prime purpose of the book is to illustrate principles and 
methods by which problems concerning the behaviour of 
systems of rigid bodies may be analysed. The calculation of 
the inertial characteristics of such systems is studied in 
Chapter 1. Principles of linear and angular momenta are 
applied in Chapter 2 to the calculation of plane motions and 
in Chapter 3 to the calculation of three dimensional motions. 
In Chapter 4, the methods associated with the names of 
Lagrange and Hamilton are used to provide alternative 
solutions of some of the earlier problems (the reader is thus 
able to compare their relative effectiveness) and are also applied 
to many new systems. The book concludes with a number of 
problems relating to oscillatory systems having two or three 
degrees of freedom. 

The topics included are those which commonly form part 
of an honours degree course in applied mathematics and it is 
hoped, therefore, that a useful supplement to the available texts 
of analytical mechanics has been provided. The exercises at 
the end of the chapters form a fairly stem test of the student's 
ability to apply the techniques illustrated in the main body of 
the book. 

In conclusion, it is a pleasure to acknowledge the assistance 
of my secretary, Mrs. Audrey Breakspear, who performed so 
efficiently the difficult task of typing the original manuscript. 

D. F. L. 
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Chapter 1 

Characteristics of Mass Distributions 


1.1 Centres of Mass If a mechanical system comprises n particles of 
masses m i (i = 1,2 situated at points P { having position vectors r, 
with respect to an origin 0, its centre of mass is the point G having position 
vector r, determined by the equation 

n n 

r = £ m,tj £ m i- ( 1 . 1 ) 

i~ 1 i = 1 

If ( x t>y l >z i ) are the coordinates of the particle m t with respect to a 
rectangular cartesian coordinate frame Oxyz and (3c, y, z ) are the coordinates 
of G, then 

x = Im, x,/L m i , etc. (1.2) 

In the limiting case when the system is regarded as a continuous 
distribution of matter over a region of space r, p(x,y,z) being the density, 
then sums are replaced by integrals, thus: 

x = | px dvj j p dv, etc., (1.3) 

r r 

the integrals being volume integrals assuming that F is a three-dimensional 
region ; if F is a surface, then p is taken to be the surface density and the 
integrals are evaluated over the surface; if r is a curve, p is the line density 
and the integrals are line integrals. 

Problem 1.1 Two distinct distributions of matter have total masses 
Af,, M 2 and their centres of mass are at the points r, , r 2 . Find the centre 
of mass of the combined distribution. 

Solution. Let m t (i = 1,2, ...,n) be the particles of the first distribution 
at the points r, and m'j(j = 1,2,..., n') be the particles of the second at the 
points tj. Then 

M i r i = Z m . r o W 2 r 2 = X m j r 'j- 

‘ i J 

If the centre of mass required is at r, then 

(M, + M 2 )7 = ]rm f r, + Zm'/j = Af,r, + M 2 r 2 . (1.4) 

This equation determines r. q 

Problem 1.2 The density inside a sphere of radius a is proportional to 
the depth beneath the surface. Determine the centre of mass ofa hemisphere. 
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Solution. Taking axes Ox, Oy in the plane face of the hemisphere, with 
0 at its centre, it is obvious that the centre of mass lies on Oz; let (0,0, z) 
be its coordinates. Introducing spherical polar coordinates ( r,9,<j>) 
(Fig. 1.1), the density is given by 

p = k(a-r). 


X 



Employing equation (1.3) and expressing the volume integrals in spherical 
polars, it follows that 

2 = jjf pz ■ r 2 sin 6 drddd<j> -r- JJJ p. r 2 sin 9 drd0d<f>. 

Since z = r cos 9, the first of these integrals reduces to the triple integral 

dr J** dQ j 2r ‘ k(a-r)r 3 sin A cos 9 d<f> = nka s /20. 

The second integral is evaluated in the same way and yields the result 
nka A !9. Thus z = -fea. □ 

Problem 1J The radii of the ends of a frustum of a uniform, solid, right 
circular cone are a and 2a. Show that the distances of its mass centre from 
the ends are in the ratio 17:11. 

X 
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Fig.1.2 


Solution. Suppose that a cone C of the same material and base radius a 
is placed on the frustum F (Fig. 1.2), thus generating a cone F+ C of base 
radius 2a. If h is the height of F, then h is also the height of C. Taking the 
x-axis as indicated, the mass centres of F, C and F + C lie on this axis and 
have x-coordinates x, 5h/4, h/2, respectively (the centre of mass of a 
uniform cone of height h is a distance from the base). If p is the density 
of the material, the masses of F and C are Ina 2 hp/3, na 2 hpfi, respectively. 
Employing equation (1.4) to give the x-coordinate of the centre of mass of 
F + C,wege, , 5 \ Ip [\ 

2 * = \ 3 X+ 3 4//(3 + 3 / 

Whence x = 1 1 /j/ 28, which is equivalent to the result stated. □ 

1.2 Moments of Inertia Let / be a straight line and let pfi = 1,2 ,...,«) 
be the perpendicular distances of particles m t (i = 1,2,.... it) from /. Then, 
the moment of inertia of the system of particles about the axis / is I, where 

I=tm iP 2 . (1.5) 

1 = 1 

If m, has coordinates (x,., y ,, z t ) in a rectangular frame Oxyz, the moments 
of inertia of the system about the coordinate axes are therefore 

A = I>,(y?+z?), B = £m t (z 2 +x 2 ), C = 'Zmfxf + y 2 ), (1.6) 

respectively. In the limiting case bf a continuous distribution of matter 
having density p(x,y,z) over a region F, 

A = jp(y 2 + z 2 )dv, etc. ( 1 . 7 ) 

r 

If M is the total mass of a matter distribution and / is its moment of 
inertia about a certain axis /, the radius of gyration of the distribution about 

/ is k, where , 

/ = Mk 2 . ( 1 . 8 ) 

Problem 1.4 Calculate the radius of gyration of the sphere described in 
Problem 1 .2 about one of its diameters. 

Solution. Taking the centre as origin and the diameter as z-axis, the 
moment of inertia is the volume integral 

fp(x 2 + y 2 ) dv, 

calculated over the sphere. Introducing polar coordinates as in Problem 
1 . 2 , x 2 +y 2 = r 2 sin 2 0 and we have accordingly to evaluate the triple 

integral f dr f" d0 [ 2 * k(a— r)r*sin 3 0 d(j>. 

Jo Jo Jo 

The result is 4ir/a 6 /45. 
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The mass of a hemisphere was found in Problem 1.1 to be tc/m*/6. 
Hence for the sphere, M = nka*/3. The radius of gyration now follows 
from equation (1.8): k = 2a/y/\5. □ 

1.3 Products of Inertia 

Problem 7.5 Calculate the moment of inertia of a matter distribution 
about an axis L, having direction cosines (/, m, n) relative to a frame Oxyz , 
and passing through 0. 

Solution. Let be a typical particle at P{x i9 y i9 z t ). PN is perpendicular 
to the axis L (Fig. 1.3). Then 

PN 2 = OP 2 —ON 2 = rf — (r i cos a) 2 . 

2 



But r. cos a = r r u = /x i + my,+ nz x , u = (/, m, n) being the unit vector in 
the direction of L. Hence the required moment of inertia is 

X m, PN 2 = 0 2 +yf+zf - ( lx,+my, + nz, ) 2 ], 

= Al 2 + Bm 2 + Cn 2 -2Fmn-2Gnl-2Hlm, (1.9) 

since l 2 + m 2 + n 2 = 1; A, B, C have been defined in the previous section 
and 

f = G = Z m f z i x i. H = (1.10) 

F, G, H are called the products of inertia of the masses m, with respect to 
the frame Oxyz. If the plane Oyz is a plane of symmetry for the matter 
distribution, then G = H = 0. If Ozx is a plane of symmetry, then 
H = F = 0 and if Oxy is a plane of symmetry, F = G = 0. □ 

For a continuous distribution, equations (1.10) are replaced by 

F = J pyz dv, etc. (1.11) 

Problem 1.6 A uniform plate is in the shape of a quadrant of a circle of 


4 


radius a. If M is its mass, calculate its moment of inertia about an axis 
through the corner inclined at an angle a to one of the straight edges. 

Solution. The moment of inertia of the plate about a straight edge is 
clearly one quarter of the moment of inertia of a complete circular plate 
about a diameter. The latter is well-known to be \M'a 2 , where M' = 4M 
is the mass of the circular plate. Thus, if we take the axes Ox, Oy along the 
straight edges, then 

A = B = \Ma 2 . 

The product of inertia with respect to the axes Oxy is 

H = | pxy dA, 

where dA is an element of area of the plate and p is its mass per unit area. 
Transforming to polar coordinates (r, 8), H can be expressed as a repeated 
integral thus: 

H = j u dr j** pr 3 sin 8 cos 8 dd 
= jpa 4 = Ma 2 /2n. 

The direction cosines of the axis are (cos a, sin a, 0). Substituting in (1.9), 
therefore, the required moment of inertia is 

^Ma 2 (cos 2 a + sin 2 a)—-Ma 2 sinacosa = Ma 2 ( 7 —^- sin 2 a \ □ 

4 n V4 2 n ) 


1.4 Theorems of Parallel and Perpendicular Axes If I G is the moment of 
inertia of a matter distribution about an axis through the centre of mass 
and / is the moment of inertia of the same matter distribution about a 
parallel axis distant d from the first axis, then 

1 = I a + Md 2 , (1.12) 

where M is the total mass of the distribution. This is the theorem of 
parallel axes for moments of inertia. 

Let GX, G Y, GZ be rectangular axes through the centre of mass of a 
distribution and let Ox, Oy, Oz be parallel axes through any point O. 
Let ( x,y,z ) be the coordinates of G with respect to the frame Oxyz. Then, 
if A g , B g , C g are the moments of inertia of the distribution about the 
axes GX, GY, GZ respectively and A, B, C are the moments of inertia 
about the axes Ox, Oy, Oz, it follows from (1.12) that 


A = A G + M{y 2 + z 2 ), 

B = B g + M(z 2 + x 2 ), - 
C = C G + M(x 2 +y 2 ). . 


( 113 ) 
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Further, if F c , G a , H G are the products of inertia relative to the frame 
GXYZ and F, G, H are the corresponding products of inertia in the frame 
Oxyz, then it may be shown that 

F = F G + M yz, G = G c + Mzx, H = H G + Mxy. (1.14) 
This is the theorem of parallel axes for products of inertia. 

If the axes Ox, Oy are taken in the plane of a thin plate (lamina) and Oz 
is perpendicular to the plate, then the moments of inertia of the plate about 
the axes are related by the equation 

C = A + B. (1.15) 

This is the theorem of perpendicular axes for a lamina. 

Problem 1.7 Ox, Oy, Oz lie along three concurrent edges of a uniform 
rectangular block of mass M. The lengths of these edges are a, b, c respec¬ 
tively. Calculate the products of inertia of the block in the frame Oxyz 
and hence find its moment of inertia about a diagonal of a face with 
dimensions a and b. 

Solution. Taking parallel axes through the centre G of the block, the 
moments of inertia about these axes are 

= + B g = ^2 M (c 2 + a 2 ), C G = ^A4(a 2 + f> 2 ). 

Since the block is symmetric with respect to these axes, the products of 
inertia all vanish in this frame. Employing the parallel axes theorems, 
the moments of inertia and products of inertia in the frame Oxyz are now 
calculated to be (3c = ja, y = $b,z = $c) 

A = \M(b 2 +c 2 ), B = $M(c 2 + a 2 ), C = $M(a 2 +b 2 ), 

F = $Mbc, G = jMca, H = \Mab. 

The diagonal of the face lying in the xy-plane has direction cosines 
(Aa,kb,0), where k = (a~ + b~) *. Substituting in equation ( 1 . 9 ), the 
moment of inertia about this diagonal is now found to be 

M(a 2 b 2 +2 b 2 c 2 +2 c 2 a 2 ) 

6(a 2 +b 2 ) ’ ^ 

1.5 Principal Axes of Inertia Let Oxyz, Ox'y'z' be two rectangular 
cartesian frames with a common origin. The transformation relating the 
coordinates of a point in the first frame to the coordinates of the same point 
in the second frame, can be written 

x' = l x , x x+l x . y y + l x . z z, j 

y' = tf*x+ly,y+lf,*, [ ( 1 . 16 ) 

z' = l z - x x+l I . y y+l t . z z, J 


6 


where l x . x is the cosine of the angle between the axes Ox' and Ox, etc. 
Denoting the column matrix with elements (x,y,z) by x and the column 
matrix with elements ( x',y',z') by x', this transformation can be expressed 
in matrix form thus: 


where 


x' = 

lx, 







1 = 


ly-y 

ly; 


U; 

K; 

l, J 


(1.17) 

(1.18) 


The elements in the first row of this matrix are the direction cosines of 
Ox' in the frame Oxyz, the second row are the direction cosines of Oy' 
and the third row are the direction cosines of Oz'. The matrix I is orthogonal 
and hence satisfies the condition 


l - 1 =I r . (1.19) 

Let A,B,C,F,G,H be the moments and products of inertia relative to the 
frame Oxyz and let A', etc. be the corresponding quantities in the ‘dashed’ 
frame. Form the matrix J thus: 


J = 


A 

-H 


-G 


-H 

B 

-F 


-G 

-F 

C. 


( 1 . 20 ) 


The matrix J' is constructed in the same way from the coefficients of inertia 
in the frame Ox'y'z'. Then 


J' = 1JI T . (1.21) 

This equation, together with the parallel axes theorems of the previous 
section, enable us to calculate the coefficients of inertia in any rectangular 
frame when these coefficients are known in one such frame. 

It is a well-known algebraic result that, given any symmetric matrix J, 
it is possible to find an orthogonal matrix 1 such that the matrix J' given 
by equation (1.21) is diagonal. After transformation to such a frame Ox*y*z*, 
the products of inertia F*, G\ H* will all vanish. The axes Ox*, Oy*, Oz* 
are then called principal axes of inertia and the moments of inertia about 
them, A*, B*, C*, are the principal moments of inertia. 

If ( l,m,n) are the direction cosines of a line through 0 with respect to 
the principal axes, the moment of inertia of the distribution about this line 
is given by 


/ = A*l 2 + B*m 2 + C*n 2 . (1.22) 

This is the special form taken by equation (1.9) for principal axes. 
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Considerations of symmetry often permit the principal axes through a 
point to be recognized without calculation. If this is not the case, we 
determine the matrix I which transforms from given axes Oxyz to principal 
axes through 0, as follows. Solve the determinantal equation 

| J—AI| = 0 (1.23) 

for A (I is the unit matrix). Let A t , A 2 , A 3 be the three roots (these are always 
real). Solve the three homogeneous equations 

(J —Ajl)x = 0 (1.24) 

for the ratios x:y:z and choose values of x,y,z to satisfy the condition 

x 2 +y 2 +z 2 = 1 . (1.25) 

These three values are l x . x> l x . y , l x . z , (i.e. the direction cosines of Ox' in 
Oxyz) and hence give the first row of 1. The second and third rows of I are 
found in the same manner after replacing A, by A 2 and A 3 respectively. 
The principal moments of inertia are given by 

A* = X t , B* = A 2 , C* = A 3 . (1.26) 

If A, = A 2 # A 3 , calculate the direction cosines of Oz' in Oxyz from the 
root A 3 and so write down the third row of 1 . The axes Ox', Oy' may be 
taken in any pair of perpendicular directions which are both perpen¬ 
dicular to Oz'. 

If A, = A 2 = A 3 , all sets of perpendicular axes through 0 are principal 
axes. 

Problem 1.8 Uniform rods OA, AB, BC, are each of mass M and length 
a and are mutually perpendicular. Determine the principal axes and 
moments of inertia at the centre of mass of this sytem. 

Solution. Taking x-axis along OA, y-axis through 0 parallel to AB and 
z-axis parallel to BC, the centres of mass of the rods have coordinates 
($a, 0,0), (a,^a,0), (a, a, {a). Let (x,y,z) be the centre of mass of the system. 
Then 

3 Mx = M({a+a+a), 3My = M({a + a), 3Mz = M .{a. 

Thus x = 5a/6, y = a/2, z = a/ 6 . 

In the frame Oxyz, the moments of inertia about the axes of the whole 
system are easily found to be 

A = fMa 2 , B = §Ma 2 , C = V Ma 2 . 

The contributions of the rods OA and AB to the product of inertia F are 
zero, since every particle of these rods has zero z-coordinate. The contribu¬ 
tion of the rod BC is 

j” paz dz = {pa 3 = {Ma 2 , 
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p being the mass per unit length. Thus 

F = {Ma 2 . 

Similar calculations yield 

G = {Ma 2 , H = {Ma 2 . 

We next apply the parallel axes theorems to give the moments of inertia 
and products of inertia with respect to a parallel frame through the centre 
of mass G. The results are 


A g = \Ma 2 , B g = {Ma 2 , C G = f Ma 2 , 
F a = {Ma 2 , G g = ±Ma 2 , H G = {Ma 2 . 


The determinantal equation (1.23) can now be constructed. It is 


~{ 


_1 

4 


12 

_JL 

4 


= 0 , 


where p = X/Ma 2 . The characteristic roots are p = i,{{ and implying 
that the principal moments of inertia at G are 


A* = Ma 2 , B* = \{Ma 2 , C* = {Ma 2 . 
Taking A = Ma 2 , equations (1.24) take the form 


x y z 
6~4~ 12 


o _f_Z_£ = o -—-l-l = o 

’ 4 2 4’ 12 4 6 


These have the solution x:y:z = 1: — 1:1. To satisfy the condition (1.25), 
we take 


x = 1A/3 y=- 1/V3 z=l/V3. 

These are the direction cosines of the principal axes about which the moment 
of inertia is A*. 

Similarly, with A = \{Ma 2 and {Ma 2 , the corresponding solutions are 
x=l/V2 y = 0, , z=—1/^2, 

X = 1A/6, y = 2/J6, z = l/ v /6. 

These are the direction cosines of the principal axes about which the 
moments of inertia are B* and C* respectively. 

The transformation matrix 1 relating coordinates with respect to the 
axes parallel to Oxyz to coordinates with respect to the principal axes is 



V 2 
V3 
_ 1 


-V 2 


0 

2 


>/ 2 ' 

-n/3 

1 . 


B 


□ 
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Problem 1.9 Calculate the principal axes and moments of inertia at a 
corner of a uniform cube. 

Solution. Taking axes through the centre of the cube perpendicular to 
the faces, the moments and products of inertia in this frame are found to be 

a g = B c = C G = %Ma 2 , F g = G g = H g = 0, 
where M is the mass of the cube and a is the length of an edge. 

Transforming to parallel axes Oxyz through a corner, we now find that 

A = B = C = f Ma 2 , F = G = H = £ Ma 2 . 

The characteristic equation (1.23) is 


_1 

4 

_1 

4 


-i 

-i 


_1 

4 

-i 

I-p 


= 0 , 


p = 


ii u 
12 ’ 12 ’ 


with (i = X/Ma 2 . The characteristic roots are 
principal moments of inertia are accordingly 

A* = ftMa 2 , B* = ftMa 2 , C* = £Ma 2 
Taking X = j£Ma 2 , the equations (1.24) all become 


and the 


x+y+z = 0 . (1.27) 

This has an infinity of solutions x:y:z. Only the third principal axis 
corresponding to the root X = %Ma 2 is unique; this root leads to equations 
(1.24) having the form 

— 2 x+y+z = 0 , x-2y+z = 0 , x-Fy-lz = 0 . 

Thus, X = y = 2 = l/y/3 are the direction cosines of this principal axis; 
it is clearly the diagonal of the cube through the corner O. Any pair of 
axes perpendicular to this diagonal can be chosen to represent the two 
remaining principal axes; all such axes have direction cosines satisfying 
equation (1.27). □ 


EXERCISES 


1. A solid is bounded by the coordinate planes and that part of the ellipsoid 

x 2 y 2 z 2 

~2+ 7T2+—J = 1 

a 2 b 2 c 2 

which lies in the first octant. If the density varies as the square of the distance 
from the origin, show that the z-coordinate of the centre of mass is given by 

_ _ 5c a 2 +b 2 + 2c 2 
16' a 2 + b 2 +c 2 ' 

(Hint: Transform by x = ar sin 9 cos (f>, y = hr sin 0 sin (f>, z = cr cost?.) 
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2. A circle of radius a is rotated about an axis / in its plane, distant r(> a) 
from its centre, to generate a torus. If the torus is constructed from uniform 
material, show that its radius of gyration about / is k, where 

k 2 =\a 2 + r 2 . 


3. A lamina lies in the xy-plane and has moments of inertia A, B and 
product of inertia H with respect to these axes. Axes Ox'y' are obtained 
from Oxy by rotation in the plane through an angle 6. Show that the 
lamina’s product of inertia with respect to the rotated axes is given by 

H’ = £(4 — B) sin 20+H cos 26. 


Hence find 9 if Ox', Oy' are principal axes. 

If the lamina is a uniform triangular plate with sides a, b enclosing an 
angle C, show that a principal axis through the vertex C makes angles 
9, <f> with the two sides such that 


tan ((/> - 9) 


(a 2 - fc 2 ) sin C 
(a 2 + b 2 ) cos C+ab' 


4. A uniform solid tetrahedron of mass M is bounded by the coordinate 
planes and the plane 


x 

a 


+S+ 


z 

c 


= 1 . 


Calculate its moments and products of inertia with respect to the frame 
Oxyz and hence show that its moment of inertia about an axis through 
O perpendicular to the opposite face is 


io M 



+ b 2 + c 2 — 


6 a 2 b 2 c 2 "I 

b 2 c 2 + c 2 a 2 + a 2 b 2 J 


5. A uniform hemisphere has mass M and radius u. Rectangular axes 
Oxyz are constructed through a point O on the circumference of its plane 
face, Oy passing through the centre of this face and Oz being tangential 
to the circular edge. Calculate the moments and products of inertia of 
the hemisphere in this frame and hence show that Oz is a principal axis 
at O and that one of the other principal axes makes an acute angle tan" ‘3 
with Ox and the other makes the same angle with Oy. Show that the princi¬ 
pal moments of inertia are ^Ma 2 , \Ma 2 , Ma 2 . 
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Chapter 2 
Plane Motions 


2.1 Systems with Few Particles If a dynamical system comprises only a 
small number of interacting particles, its motion can often be calculated by 
writing down the equations of motion of the particles separately. These 
equations of motion follow from Newton’s second law, which is only valid 
for accelerations calculated relative to an inertial frame, i.e. a set of axes 
whose origin is moving uniformly relative to the centre of our galaxy 
and which is non-rotating relative to the galaxy. Unless otherwise stated, 
all dynamical quantities (e.g. acceleration, linear momentum, etc.) which 
arise in the problems which follow, will be supposed calculated within an 
inertial frame. However, a frame fixed in the earth will often be treated as 
if it were inertial, with the consequence that the results calculated in such 
cases will be approximate only. 

Problem 2.1 A particle P of mass m lies on a smooth horizontal table 
and is attached to a long, inextensible string which passes through a 
smooth hole in the table. The other end of the string carries a particle Q of 
mass km hanging freely. The particle P is projected from rest along the table 
at right angles to the string with speed y/(%ag), when it is a distance a 
from the hole. Prove that the particle Q will begin to rise if k < 8 , and 
that, if the greatest distance of P from the hole is 2a, then k = 3. 



Solution. Let (r, 0) be the polar coordinates of P at time t with respect to 
the hole as pole (Fig. 2.1). Then, if T is the tension in the string, resolving 
the equation of motion for P in the radial and transverse directions, we 

° b,ai ” rtf-rO") = -T, (2.1) 


- 0 . 

r dt 


( 2 . 2 ) 
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Since the string is inextensible, the acceleration of Q is r vertically upwards 
and its equation of motion is 

kmr = T—kmg. (2.3) 

Integrating equation (2.2) and employing the initial conditions, we get 
r 2 6 = constant = aj(iag). (2.4) 

Eliminating T and & from equation (2.1) by using equations (2.3), (2.4), 
the following equation is obtained: 

(k+l)r = (Za 3 g/r 3 )-kg. (2.5) 

Initially r = a and it follows that at this instant r = (%-k)g/(k+ 1 ). 
Thus r > 0 initially, and Q commences to rise provided k < 8. 

Since Y = d(\r 2 )/dr , equation (2.5) integrates immediately to yield 

j(k + l)r 2 = (k+A)ag - kgr - (4a 3 g/r 2 ), 

the constant of integration following from the initial conditions r = a, 
r = 0. If P is at its maximum distance from the hole when r = 2a, r must 
vanish for this value of r also. Thus k = 3. □ 


Problem 2.2 Let m lt m 2 be the masses of a planet and its satellite (or 
the masses of the components of a binary star system), attracting one 
another according to the Newtonian law of gravitation. Discuss their 
motions. 

Solution. If r,, r 2 , are the position vectors of the centres of these two 
bodies with respect to the origin 0 of an inertial frame, the forces F 1; F 2 
exerted upon them are given by 


Fi = 


Gm 1 m 2 (r 2 -r 1 ) 


P 2 = 


Gm l m 2 (r) —r 2 ) 


l*l -, 2l 1*2 *21 

where G is the gravitational constant (Fig. 2.2). Evidently F, = — F 2 by 
Newton’s third law. The equations of motion of the bodies are therefore 


Gm 1 m 2 (i 2 -T l ) 

ki-^l 3 ’ 


m 2 *2 = 


i 13 • 



Fig. 2.2 
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Adding these equations, we find that 

m \ r i + W 2 ^2 = 0 

and hence, if r is the position vector of the centre of mass of the bodies, 
it follows that 

(m l +m 2 ) r = 0. 

Thus, the centre of mass has zero acceleration and moves with uniform 
velocity in any inertial frame. 

Let r = r 2 — r, be the position vector of the satellite in a frame S whose 
axes are parallel to the original frame and whose origin is at the centre 
of the planet. It then follows from equations (2.6) that 


f = f 2 -f, = -G(m, + m 2 )r/r 3 , (2.7) 

where r = |r|. This is exactly the equation of satellite motion which would 
be derived if S were treated as inertial (i.e. the planet were taken to be 
unaccelerated) and the mass of the planet were taken to be (mj+m 2 ); 
this we shall term the corrected mass of the planet. 

Taking the vector product of both sides of equation (2.7) with r, it 
follows that 

d(r x r)/dt = r x f = 0. 

Thus r x f = h, (2.8) 

where h is a constant vector, and this implies that r is always perpendicular 
to h, i.e. that the satellite’s trajectory lies in a fixed plane through the 
planet. 

Taking polar coordinates ( r,6 ) in the plane of the satellite’s motion, the 
planet lying at the pole, the radial and transverse components of equation 
(2.7) are (fi = G(m, + m 2 )) 

r- r p = il (r 20) = 0. (2.9) 

Thus r 2 6 = h (a constant), (2.10) 

which is a scalar form of equation ( 2 . 8 ). h is the moment of the satellite’s 
velocity about the pole. Putting u = 1/r, then 0 = hu 2 and hence 


. _ 1 du 1 du . du 

r ~ ~^~dt ~ ~ ~ 7a' 


u 2 d6 


dO ’ 


d ( du\dO , 2 2 d 2 u 
r = —— h— — = -nu— 
d0\ dOJ dt dd 2 


The first of equations (2.9) can now be written in the form 

d 2 u fi 


This is a linear equation with constant coefficients whose general solution 
is 

u = A cos (6 — a) + (n/h 2 ), 

A, a being constants of integration. By appropriate choice of the line 
6 = 0, a can always be reduced to zero. Thus this equation is equivalent to 

l/r = 1 + e cos 0 , ( 2 . 11 ) 

where / = h 2 /n and e = Ah 2 foi. Equation (2.11) is the polar equation of a 
conic of eccentricity e and semilatus rectum / with a focus at the pole. 
It has accordingly been shown that the orbit of the satellite as observed 
from the planet is an ellipse (e must be less than one, for otherwise the 
satellite recedes to an infinite distance along a parabolic or hyperbolic 
path), with a focus at the planet’s centre. 

Integrating equation (2.10) over a time interval T during which the 
satellite completes one circuit of its orbit from 6 = 0 to 6 = 2n we obtain 

2A = If fl d0 = hT ’ 

where A is the area of the ellipse. But, A = nab , where a, b, are the semi¬ 
axes of the ellipse and h = (/i/)* = i u*b/a*. It follows that 

T = 2 na*/n*. 

This equation permits a calculation of from observed values of the periodic 
time T and the semi-major axis a . It also leads to Kepler’s third law, 
T 2 x a 3 . □ 

2.2 Conservative Systems. Energy Equation The principle of work for a 
particle states that the net work done by the forces acting upon it is equal 
to the increase in its kinetic energy. This principle can be extended to any 
system of particles, yielding the result that the net work done by the forces 
of the system (both internal and external) is equal to the increment in the 
overall kinetic energy of the system. 

In certain cases, the work done by the internal and external forces as the 
system moves from one configuration to another, depends only upon the 
end configurations. In these circumstances, the potential energy of the 
system in a given configuration is defined as the work which would be done 
by the forces of the system if it were to move from the given configuration 
to some standard configuration (which can be chosen arbitrarily). For 
such a system, the principle of work can be expressed as the principle of 
conservation of energy , 

kinetic energy 4 * potential energy = constant. ( 2 . 12 ) 
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The kinetic energy of a moving rigid body of mass M can be calculated 
in two parts, (i) the kinetic energy of the overall translatory motion, which 
is that of a particle of mass M moving with the centre of mass of the body, 
and (ii) the kinetic energy of the rotary motion relative to axes through the 
centre of mass parallel to the axes of the inertial frame being employed. 
If the motion of the rigid body is parallel to a plane n and v is the velocity 
of its centre of mass, w is its angular velocity about an axis through the 
centre of mass perpendicular to n, I G is its moment of inertia about this 
axis, then its kinetic energy is given by the formula 

T = iMv 2 +jI G co 2 . (2.13) 

In the special case of a rigid body rotating about a fixed axis, its kinetic 
energy can be calculated immediately from the formula 

T = (2.14) 

where / is the moment of inertia of the body about the axis and to is its 
angular velocity. 

Problem 23 A uniform rod is placed with one end on a smooth vertical 
wall and the other on a smooth horizontal ground. It is released and slides 
in a vertical plane. Determine its motion. 



Solution. Let the smooth reactions of the wall and ground on the rod be 
P and Q as shown in Fig. 2.3. The weight Mg acts through G the centre of 
the rod and is the only force which does work as the rod moves. The weight 
is a conservative force and the work done by this force as the rod moves from 
its position in the figure to a horizontal position on the ground (the stan¬ 
dard configuration) is Mga cos 9 (2 a is the length of the rod). This, therefore, 
is the potential energy of the rod. 

A semicircle can be constructed on AB as diameter and will pass through 
0; hence OG = GA = GB = a. Thus, the triangle AGO is isosceles and 
angle GOA is 9. It follows that G describes a circle, centre 0, radius a, 
with angular velocity 0. The angular velocity of the rod about G is also 0 
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and its moment of inertia about this point is $Ma 2 . The kinetic energy of 
the rod is now calculated from equation (2.13) to be 

%M(ad) 2 +\.%Ma 2 .0 2 = \Ma 2 6 2 . 

The energy equation for the rod’s motion is therefore 
§Ma 2 0 2 + Mga cos 9 = constant. 

Assuming that the rod’s initial inclination to the wall is a, we have 0 = 0 
when 9 = x. This fixes the constant in the last equation, and hence 

0 2 = 3g(cos a — cos 9)/2a. (2.15) 

This differential equation determines the motion so long as the rod remains 
in contact with the wall (see Problem 2.5). It can be integrated to yield 9 
as a function of t, but the process leads to an elliptic integral. □ 

Problem 2.4 A uniform circular cylinder of mass M and radius a lies 
at rest on a perfectly rough horizontal plane. A particle of mass m is fixed 
to it at a point on its highest generator. If the equilibrium is slightly dis¬ 
turbed, calculate the subsequent motion. 




Fig. 2.4 


Solution. The external forces acting upon the system (cylinder + particle) 
are (i) the weight Mg of the cylinder acting at the centre of mass 0, (ii) 
the weight mg of the particle P, (iii) the normal component R of the reaction 
of the plane and (iv) the frictional component F of this reaction. Of these, 
only the weight of the particle does work during the motion; F does no 
work, provided the cylinder does not slip, since the cylinder rotates about 
C as instantaneous centre until a new particle of the cylinder is brought 
into contact with the plane and acts as instantaneous centre in its turn; thus 
the point of application of the frictional force is always stationary. Since 
weight forces are conservative, the potential energy of the system in the 
configuration shown in Fig. 2.4 is mga cos 9. The internal forces, of course, 
do no work. 

When the cylinder has rotated through an angle 9, its angular velocity 
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about 0 will be 0 and the velocity of O will be ad horizontally. Thus, the 
kinetic energy of the cylinder is 

$M(a6) 2 +l±Ma 2 d 2 = f Ma 2 0 2 , 

\Ma 2 being its moment of inertia about 0. The particle has the velocities 
a6 perpendicular to OP and ad horizontally, simultaneously; the triangle 
of velocities (Fig. 2.4) shows that its resultant velocity is v, where 

v 2 = 2 a 2 0 2 (l + cos 0 ). 

The energy equation for the system can now be written down: 

%Ma 2 d 2 + ma 2 d 2 (l + cos 9)+mga cos 6 = const. 

Initially 0 = 0, d = 0 and hence const. = mga. The equation of motion 
is accordingly 

{|M+m( 1 + cos 6)}ad 2 = mg{ 1 — cos 0 ). (2.16) 

A further integration introduces elliptic integrals. □ 

2.3 Motion of the Centre of Mass The centre of mass of a system of 
particles moves as if the whole mass of the system were concentrated there 
and ail the external forces were applied there. Thus, if M is the total mass 
of the system and F £ are the external forces, then 

£ F £ = Mr = M dv/dt. (2.17) 

The motion of the centre of mass is completely independent of the internal 
forces between particles of the system. 

Problem 2.5 Calculate the forces exerted by the wall and ground on the 
rod in Problem 2.3. 

Solution. The horizontal distance of G from the wall (Fig. 2.3) is a sin 0. 
The horizontal component of the acceleration of G is accordingly 

d 2 (a sin 0)/dt 2 = a(cos 00 -sin 00 2 ). 

Resolving the external forces acting on the rod horizontally, therefore, 
we get 

P = Ma(cos00-sin00 2 ). (2.18) 

Differentiating equation (2.15) with respect to t, we find that 
200= 3g sin 00/2a. 

We can now substitute for 0 and 0 2 in equation (2.18) and obtain the result 
P = iMg(3cos0—2costx)sin0. 

This equation shows that the rod loses contact with the wall when 
c °s0 = | cos«, i.e. when the end A has descended to two thirds of its 
initial height above the ground. 
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r The vertical acceleration of G (upwards) is 

d 2 (a cos 8)/dt 2 = — a(sin 00 + cos 00 2 ) 

= “|^(1+2 cos a cos 0-3 cos 2 0). 
Resolving the external forces vertically, this leads to the result 
Q = i Mg( 1 - 6 cos a cos 0 + 9 cos 2 0 ). 

It is easy to verify that Q is always positive and that the instant the rod 
leaves the wall, Q = \Mg. □ 

2.4 Conservation of Momentum In the case where the vector sum of the 
external forces acting upon a system has zero component in a certain 
direction, if this direction is taken to be that of the x-axis, equation (2.17) 
yields 


Mx = 0 . 

(2.19) 

Integrating, it follows that 


Mx = constant. 

( 2 . 20 ) 

Equation (1.2) implies that this is equivalent to 


p x = Mx = X m i *1 — constant. 

( 2 . 21 ) 


where p x is the x-component of the net linear momentum p of the system. 
This equation is referred to as the principle of conservation of linear 
momentum. 


Problem 2.6 Two equal uniform rods AB , BC , each of mass m, are 
smoothly jointed together at B and lie along a straight line A BC on a smooth 
horizontal plane. A particle of mass 4 m is connected to B by a light string 
which is just taut and lies along a straight line perpendicular to ABC 
also lying in the plane. The pivot B is held stationary and the rods are given 
angular velocities a) in opposite senses, so that their ends A, C commence to 
move away from the particle 4m. The pivot is released and the system 
moves on the plane. Calculate the tension in the string when the rods have 
rotated through an angle 0 . 


A 
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Solution. The motion is evidently symmetrical about the line of the 
string. In the configuration shown in Fig. 2.5, let v be the velocity of the 
pivot B and, hence, of the mass 4m. Then the centre of mass of each rod 
has velocity ad relative to B ; the velocity of the centre of mass relative to 
the plane is accordingly found by summing this velocity vectorially with 
the velocity of B as indicated. 

Since no external forces act on the system in its plane of motion, its 
linear momentum is conserved. The component of the linear momentum 
of a rod parallel to the string is found from equation ( 2 . 21 ) to be m(v + 
a6 cos 0 ) (i.e. mass of rod times the component of the velocity of the centre 
of mass parallel to the string). We can now write down an equation expres¬ 
sing conservation of momentum for the system in this direction: 

4mv + 2m(u + ad cos 0) = constant. 

Since, initially, v = 0, 9 = 0, d = co, this gives 

v = %a((o — d cos 0). ( 2 . 22 ) 

The kinetic energy of a particle of mass m moving with the centre of mass 
of a rod is 

\m(v 2 4- a 2 d 2 + 2 adv cos 0 ). 

The kinetic energy of the rotational motion of a rod relative to its centre 
of mass is %ma 2 d 2 . We can now construct the energy equation for the system 
in the form 

^.4mt ; 2 + m(tf 2 + a 2 d 2 + 2adv cos 9)+$ma 2 d 2 = const. 

Utilizing the initial conditions, this gives 

3v 2 + %a 2 d 2 + 2advcos9 = fa 2 co 2 . (2.23) 

Equations (2.22), (2.23), completely fix the behaviour of the system. 
Eliminating v between them, we find that 

d 2 = co 2 /(l+$sin 2 0). (2.24) 

Then, substituting in equation (2.22), an equation for v emerges, 

v = ^acu[l - cos 0/(1 + ^ sin 2 0)*]. (2.25) 

Differentiating equation (2.25), we get 

v = 4 aco 2 sin 0/(3 + sin 2 0) 2 . 

But if T is the tension in the string, the equation of motion of the particle 4m 
is 7 = =4 mv. Thus, a formula for Tin terms of 0 now follows immediately. 

□ 

Problem 2.7 A ring of mass m is attached to one end of a uniform rod 
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r of length 2a and mass 2m and is then threaded on to a smooth rigid 
horizontal wire. The rod is set swinging in the vertical plane containing 
the wire. Discuss the motion. 

Solution. The external forces acting upon the system (rod + ring) are 
the two weight forces mg, 2mg and the smooth reaction R of the wire on 
the ring. Since all these forces are vertical, the horizontal momentum of the 



system will be conserved. Taking the velocity of the ring to be v and the 
angular velocity of the rod to be d (Fig. 2.6), the equation expressing this 
principle is 

mv -f 2m{v - ad cos 0) = constant. 

We shall suppose that, initially, 0 = a, d = 0 and v = 0. Then, this equa¬ 
tion gives 

v = § ad cos 0. (2.26) 

Conservation of energy for the system is expressed by the equation 
%mv 2 + m(v 2 + a 2 d 2 — 2adv cos 9)+\ma 2 d 2 — 2mgacosQ = constant. 
Making use of the initial conditions, this reduces to 

9v 2 + 8 a 2 d 2 — 1 2adv cos 0 = 12^a(cos 0 — cos a). (2.27) 

Equations (2.26), (2.27) determine the motion. Eliminating v , we get 

a0 2 (l + sin 2 0) = 30(cos0-cosa). (2.28) 

Thus, if a and 0 are small, the oscillations of the rod are governed by the 
approximate equation 

d 2 = 3g(ct 2 — 9 2 )/2a, (2.29) 

all but second order terms having been neglected. Equation (2.29) indicates 
that small oscillations are synchronous with similar oscillations of a 
simple pendulum of length 2a/3. 

To calculate R , it is necessary to resolve equation (2.17) for the system 
vertically. In this case, it is convenient to calculate the right-hand member 















in two parts, by application of equation (1.4); thus, differentiating this 
equation twice, we obtain 

Mf = M l i l + M 2 r- 2 , (2.30) 

where M,, Af 2 denote the masses of the ring and rod respectively and 
r, , f 2 are the respective accelerations of the ring’s and rod’s centres of 
mass. Since the ring has no vertical acceleration, the equation for R takes 
the form 

R-3mg = 2m(ad 2 cos d + ad sin 0). (2.31) 

By differentiating equation (2.28), d can be calculated thus 

ad = 3g sin 0(2 cos 6 cos a+sin 2 0 — 3)/2( 1 + sin 2 0) 2 . (2.32) 

Then, substituting from equations (2.28), (2.32) into equation (2.31), it is 
found that 

R = 3mg(cos 2 6 - 2 cos a cos 6+ 2)/( 1 + sin 2 0 ) 2 . (2.33) □ 

2.5 Equation of Angular Momentum The sum of the moments of the 
linear momenta of a system of particles about a point C, is called the 
angular momentum h of the system about C. If r, is the position vector of 
the particle m, relative to C and v, is the velocity of this particle in the 
inertial frame F being used, then 

h = £r,xro ( v,. (2.34) 

If C is the origin 0 of F, then f, = v, and it follows that 
— = x r, x m,v ( + X r, x m j v i 

= £ T i * (F« + Fi/) = X r, x F*, (2.35) 

where F fE , F l7 are the respective external and internal forces acting on m t . 
Thus, the rate of change of the angular momentum about 0 is equal to the 
sum of the moments of the external forces about O. 

Like the energy, h can always be expressed as the sum of two parts, 
(i) the angular momentum about 0 of a particle moving with the centre 
of mass and of mass equal to that of the whole system and (ii) the angular 
momentum about the centre of mass of the motion relative to axes through 
the centre of mass parallel to the axes of F. Thus, we write 

h = h c + rxMv, (2.36) 

where r is the position vector of the centre of mass relative to O, v is the 
velocity of the centre of mass relative to F and h G is the part (ii). If the 
system is a rigid body, its motion relative to the centre of mass is a rotation 
with some angular velocity w. 
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In the special case of plane motion, h is always normal to the plane of 
motion n and only its magnitude need be considered. For a rigid body in 
plane motion, the angular velocity of its motion relative to the centre of 
mass is also normal to n and 


h G = I G (0, (2.37) 

where I G is the moment of inertia of the body about the axis normal to n 
through G. 

Differentiating equation (2.36) and noting that t = v, it is found that 
h = h c + rxMv = h c + Xfx F,e- It now follows from equation (2.35) 
that 

1 »g = X ( r i - f) x f .e = M c , (2.38) 

where M G is the sum of the moments of the external forces about G. For 
the special case of plane motion of a rigid body, equation ( 2 . 37 ) shows that 
equation (2.38) takes the form 

I G (b = M c . (2.39) 

Problem 2.8 A uniform cylinder of mass m and radius a rolls inside a 
fixed cylinder of radius b with its axis horizontal; the cylinders are in 
contact along a common generator. Discuss the motion. 



Solution. At time t, let the plane containing the generator of contact and 
the axis of the moving cylinder be inclined at an angle d to the vertical 
(Fig. 2.7). Then, if G is the centre of mass of the moving cylinder, its 
velocity is (b-a)0 perpendicular to CG. If w is the angular velocity of this 
cylinder, the velocity of a point of contact C relative to G is aw perpendicular 
to CG. Thus the velocity of C relative to a fixed observer is (b-a)6-aw. 
This must vanish and hence 

aw = ( b-a)0. 

Let (F, R) be the frictional and normal components of the reaction at C. 
Taking moments about G and using equation (2.39), we obtain 

— Fa = jma 2 w. 

Hence F = -$m(b-a)'d. 
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The acceleration of G has a component ( b-a)6 2 along CG and a com¬ 
ponent (b — a)6 perpendicular to CG. It follows that 

R-mg cos 6 = m(b — a)9 2 , (2.40) 

F—mg sin 6 = m(b—a)9. 

Hence F = %mg sin 0, 

d = - 2g sin 9/3(b — a). 

The last equation shows that the cylinder can oscillate about its lowest 
position in perfect synchrony with a simple pendulum of length 3(b-a)/2. 
Since 0 = d(j0 2 )/dO, we can integrate the last equation to give 

6 2 = fl 2 - 4g( 1 — cos 9)/3{b —a), 

where we have supposed 9 = Q, 6 = 0, initially. It now follows from 
equation (2.40) that 

R = ^mg(7 cos 0 — 4)+ m(b—a)Q 2 . 

If the moving cylinder is to make a complete circuit of the fixed cylinder, 
both 9 2 and R must be positive for 9 = n. The last two equations indicate 
that the condition for this is fi 2 > 1 lg/3(b — a). □ 

Problem 2.9 A uniform circular cylinder is placed on a horizontal plane 
and an identical cylinder is balanced on top of it so that the axes of the 
cylinders are parallel. The equilibrium is slightly disturbed. Determine the 
motion if there is no slipping. 
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Solution. At time t, let 4> be the angle made by the plane containing the 
axes of the cylinders with the vertical and let 9 be the angle through which 
the lower cylinder has rotated. Suppose X, Y are points, one on each 
cylinder, which were in contact in the equilibrium configuration. Then, 
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if C is the point of contact at time t, 

angle YAC = angle XOC = (f> —6. 

It follows that the angle made by A Y with the downward vertical is2<f> — 9 
and that the angular velocity of the upper cylinder is 2<j> - 9. 

If a is the radius of either cylinder, relative to 0, A moves in a circle of 
radius 2a with angular velocity <p. O ’s velocity is ad horizontally. The 
velocity v of A in a fixed frame is obtained by compounding these two 
velocities as indicated in Fig. 2.8; thus 

v 2 = a 2 (9 2 + 4<f> 2 +40(j> cos <f>). 

The energy equation of the system can now be constructed thus: 

$ma 2 9 2 + ±ma 2 9 2 + jma 2 (9 2 + 4(j> 2 + 4 6<j> cos <f>) 

+ {ma 2 (2(f> — d) 2 + 2mga cos <f> = constant. 

Since, initially, <j> = 0, 6 = 0 = 0, this equation reduces to 

3d 2 + 6<f > 2 + 29<p{2 cos (t> — 1 ) = 4g( 1 — cos <j>)/a. (2.41) 

Let P be the frictional component of the interaction between the cylinders 
at C. Taking moments about A for the upper cylinder, we obtain the 
angular momentum equation 

Pa = jma 2 (2$-0). (2.42) 

If F is the frictional component of the reaction of the plane upon the 
lower cylinder, by taking moments about O, an angular momentum 
equation for the lower cylinder is found: 

Pa—Fa = jma 2 9. (2.43) 

The centre of mass of the two cylinders is at C and, relative to O, this 
point moves on a circle of radius a and centre 0. It follows that C has 
accelerations (i) a<j> 2 along CO, (ii) a<}> perpendicular to OC and (iii) a'9 
horizontally. Resolving the external forces acting upon the whole system 
horizontally and employing the ljnear momentum equation (2.17), we 
now get 

F = 2m(a9+a<j> cos <l>—a<j> 2 sm<t>). (2.44) 

Eliminating P and F between equations (242H2.44), the following 
equation is derived, 

30 = <£(1 —2cos </>) + 2<i> 2 situt> = ^[0(1-2 cos </>)]. (2.45) 

Integration yields the equation 

3d = 0(1 — 2 cos <t>). (2.46) 

25 




c 












t 




So long as the cylinders remain in contact, the system’s behaviour is 
governed by the equations (2.41), (2.46). By eliminating 6 between these 
equations, a first order equation for tp is obtained, 

0 2 (17+4cos$-4cos 2 4>) = 12^(1 —cos <f>)/a. 

By resolving horizontally for the external forces acting upon the lower 
cylinder, an equation for the normal reaction Q between the cylinders at 
C can also be found and, hence, the value of <P when the cylinders separate. 

□ 

2.6 Conservation of Angular Momentum In the event that the external 
forces acting upon a dynamical system have zero net moment about a 
fixed centre C, equation (2.35) yields the integral 

h = constant. (2.47) 

This is the principle of conservation of angular momentum. 

If the external forces have no net moment about a certain axis, e.g. 
the z-axis, the component of angular momentum about this axis is 
conserved: 

h z = constant. (2.48) 

Problem 2.10 A light smooth tube can rotate freely in a horizontal plane 
about its mid-point 0 which is fixed. A uniform rod of mass m and length 
2a, is placed inside the tube so that its mid-point M is a distance a from 0. 
The system is given an initial angular velocity co about the vertical through 
0. Discuss the motion. 

Solution. At time t, let 0 be the angle through which the tube has rotated 
and r the distance of M from O. Since the external forces acting on the 
system, tube + rod, have no moment about a vertical axis through 0, 
the angular momentum of the system about this axis is conserved. The 
angular momentum of the rod about O can be calculated from equation 
(2.36) in the form jma 2 () + mr 2 6. 
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The angular momentum of the tube is negligible, so that the equation 
expressing conservation of angular momentum is 

(a 2 + 3r 2 )& = 4 a 2 oj. 


The principle of conservation of energy leads to the equation 
jm(r 2 + r 2 6 2 )+%ma 2 6 2 = § ma 2 co 2 . 

These two equations determine the motion. Eliminating 0, we find that 


r = law 



indicating that ultimately, as r -► oo, r -> law/Jl, i.e. the rod moves 
uniformly relative to the tube. n 


EXERCISES 

1. Two particles P and Q, each of mass m, are connected by a light elastic 
string of modulus X and natural length a. P is placed close to the edge of a 
smooth horizontal table and Q is placed on the table so that the string is 
just taut and perpendicular to an edge. P is pushed gently over the edge. 
If P, Q move distances x, y respectively in time t, show that 

x +y = 2 9t 2 , z + n 2 z = g, 

where z = x— yandn 2 = 2X1 ma. Show that, when the string first becomes 
slack, Q will still be on the table provided 2X > n 2 mg. 

2. A block of mass m is placed with a plane face in contact with a smooth 
horizontal plane. A circular cylindrical cavity of radius b, with a horizontal 
axis, exists within the block and a uniform solid circular cylinder of 
radius a and mass m is free to roll without slipping on the wall of the cavity 
with its axis parallel to that of the cavity. The system is set in motion so that 
the solid cylinder executes complete revolutions of the cavity and the block’s 
motion is perpendicular to the cavity axis. If the net linear momentum of 
the system is zero and contact between the cylinder and the cavity is about 
to be broken when the cylinder is at its highest point, show that 

6 2 (b—a)(2 + sin 2 0 ) = 2^(3 + 2 cos 0), 

where 0 is the angle made by the plane containing the axes of the solid 
cylinder and cavity with the downward vertical. Show that the normal 
and frictional components of the reaction between the cylinder and the 
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cavity wall are given by 

3mg{\ + cos 0X3 + cos 0) _ /M 0 sin 0(3 + 3cos0+cos 2 0) 

R= ( 2 +sin 2 0) 2 ’ ” (2 + sin 2 0) 2 * 

3. Two uniform rods AB and BC, each of length 2 a, and having masses 
m, 2m respectively, are smoothly hinged together at B. The rods can move 
in a vertical plane, the ends A and C sliding on a smooth horizontal plane. 
If 0 is the inclination of AB to the horizontal and 0 = a initially, obtain 
the equation of motion 

a6 2 (\2 — sin 2 0 ) = 18g(sina — sin 0 ). 

4. A uniform circular cylinder can turn freely about its axis, which is 
fixed horizontally. A second identical cylinder is placed with a generator 
lying along the highest generator of the first cylinder, and the system is 
given a slight disturbance from rest. If the coefficient of friction between the 
cylinders is n, show that slipping begins when the angle 0 between the 
plane of the axes and the upward vertical is such that 

13/icos0—sin0 = 8 /i. 

5. A bead P of mass m is threaded on a smooth, uniform circular wire of 
mass m and radius a, free to rotate in a horizontal plane about a smooth 
vertical axis through a point 0 of its circumference. With the wire station¬ 
ary, the bead is projected with velocity u from the point A, where OA is a 
diameter of the wire. If 0 is the angle PC A, where C is the centre of the 
wire, obtain the equation of motion 

a 2 0 2 (3-cos 2 0) = 2m 2 cos0. 
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Chapter 3 

Three Dimensional Motions 


3.1 Euler’s Equations Consider a rigid body which rotates about a fixed 
point 0 of itself as pivot. Choosing a rectangular cartesian reference frame 
Oxyz, let to = (co x ,a) y ,co z ) be the angular velocity of the body about 0 
and let h = {h x ,h y ,h z ) be its angular momentum about 0. Then, if J is 
the inertia matrix (equation 1.20) for the body with respect to the frame 
Oxyz, it can be proved that 


h = Jto, (3.1) 

or h x = Ato x — Hto y —Gto z , 

h y = - Hto x +Bw y — Fto z , 
h t = -G(o x —F<o y + Cto z , 

where the components of h and to are to be arranged as columns. In the 
special case when the axes of the frame are principal axes, F = G = H = 0 
and equation (3.1) gives 

K = Aa> x> h y = Bo) y , h t = Cto z , (3.2) 

where A, B, C are now the principal moments of inertia. 

Suppose that the frame Oxyz moves with the body, so that its axes are 
always principal axes. Then the frame has angular velocity to relative to 
the inertial frame instantaneously coincident with it. Thus, as proved in 
section 6.4 of Vector Algebra (Marder) in this series, 

dh oh 

-r = —+<oxh, (3.3) 


where dh/dt denotes the rate of change of h relative to the inertial frame and 
dh/dt denotes the rate of change of h relative to the rotating frame. Let M 
be the moment about 0 of the external forces applied to the body. Then 
equations (2.35), (3.3) imply that 

{dh/dt)+to xh = M. (3.4) 


Substituting for h from equations (3.2), we derive Euler's equations, 
A(b x + {C—B)a> y (o z = M x , ' 
Btb y +{A-C)to z (o x = M y , ■ 

C(b z +{B — A)a) x (o y = M z . 

Problem 3.1 Discuss the motion of the earth about its centre. 


(3.5) 
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Solution. The earth’s centre can be treated as a stationary pivot O. Taking 
Oz through the poles and Ox, Oy through points on the equator, since the 
earth is very nearly an oblate spheroid, these will be principal axes with 
A = B and C > A. Neglecting tidal forces due to the sun and moon, the 
motion about 0 takes place under zero external forces. Thus Euler’s 
equations give 

Aw x +(C—A)w y w z = 0, Aw y — (C—A)w z w x = 0, Cw z = 0. 

Thus w z = n (a constant) and 

w x + vw y = (b y — vo) x = 0 

where v = n(C — A)/A. These equations have the general solution 
co x = a cosfvt + a), a) = asin(vt + a). 

Thus 

h x = /4acos(vt + a), h y = Aa sin(vf + a), h z = Cn. 

These equations show that the angular momentum vector h can be repre¬ 
sented by OP in Fig. 3.1; in this diagram, ON = Aa, PN = Cn and ON 

2 



makes an angle vt + a with Ox. It is now clear that OP traces out a cone with 
Oz as axis and semi-vertical angle <5, where tan 5 = Aa/Cn. But, since the 
motion takes place under no external forces, equation ( 2 . 35 ) shows that h 
is a constant vector in a non-rotating frame, i.e. OP is a fixed direction in 
space. Relative to a non-rotating frame, therefore, the earth’s axis Oz 
rotates about a fixed axis OP, thus generating a cone in space with semi¬ 
vertical angle S. This is the phenomenon of Eulerian nutation. The rate of 
rotation is evidently v, so that a complete cycle is completed in time 

2n/v = 2nA/n(C — A). 

As the earth’s axis describes its cone, the earth rotates about this axis 
with angular velocity n. □ 
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Problem 3.2 Discuss the motion of a body moving freely about a fixed 
pivot 0 under no other forces. 

Solution. Multiplying Euler’s equations (3.5), with M x = M y = M z = 0, 
by co x ,o) y , oj z respectively and adding, we get 

Aco x (b x + Ba) y (i} y +C(o z (b z = 0. 

Upon integration, this yields 

Awl + Boi y + Coj? = 2T, (3.6) 

where T is a constant which will later be identified with the kinetic energy. 

Again, multiplying Euler’s equations by Aw x , Bw y , Cw z , respectively 
and adding, we obtain 

A 2 a> x (b x + B 2 w y (b y + C 2 (o z (b z = 0. 

This integrates to give 

A 2 wl + B 2 (ol + C 2 col = H 2 . (3.7) 

H is the magnitude of the angular momentum vector h (see equations 3.2). 

Equations (3.6), (3.7), can now be solved for (o x ,w z in terms o(w y ; then, 
substituting in the second Euler equation, a first order differential equation 
for w y is obtained. In general, the solution of this equation involves elliptic 
functions, however, if the initial conditions are chosen appropriately, a 
solution in elementary terms can arise. For example, if H 2 = 2 BT, it will 
be found that 


Awl _ 2 T-Bwl _ C( °l 
C-B C-A ~B-A ’ 


(3.8) 


where we are assuming A < B < C. These equations permit us to reduce 
the second Euler equation to the form 

w y = a(A 2 -a$, 

where a = \_(B - A\C - B)/AC]*, A 2 = 2 T/B. 

Choosing the origin of time appropriately, this equation integrates to 
yield 


w y = AtanhAat. 

Substituting in equations (3.8), this leads to 


AsechAaf> - [H^]^ sechAaf - 

Clearly, as t -» cc,w y -► A and w x -» 0 , w y —► 0, i.e. ultimately, the body 
rotates about the principal axis of inertia corresponding to the inter¬ 
mediate moment of inertia B, with uniform angular velocity A. □ 
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3.2 Kinetic Energy of a Rigid Body If the body rotates about a fixed 
pivot 0 relative to an inertial frame Oxyz with angular velocity to, its 
kinetic energy will be given by the equation 

T = i/co 2 , (3.9) 

where I is the body’s moment of inertia about the instantaneous axis of 
rotation. If this axis has direction cosines ( l,m,n ), / is given by equation 
(1.9). Since to x = Ito, to y = mto, co z = ncu, equation (3.9) can be written 
in the form 

T = %Ato 2 + Bw 2 + Cto 2 - 2Fco y to z - 2 Geo, to x - 2 Hco x (o y ). (3.10) 

In the special case when the axes Ox, Oy, Oz, are principal axes, this 
formula reduces to 

T = %Aco 2 x + Bto 2 + Cto 2 ). (3.11) 

The kinetic energy of a rigid body moving generally can now be cal¬ 
culated by separating this into two parts, (i) the kinetic energy of the 
translatory motion, and (ii) the kinetic energy of the rotary motion about 
the centre of mass, as explained in section 2.2. 

Problem 3J A uniform square plate has one comer freely pivoted to a 
fixed point 0 of a smooth horizontal plane. The plate is initially vertical 
and then falls from rest under gravity, its lower edge sliding on the plane. 
Calculate the motion. 





Solution. Let Ox be the initial position on the plane of the edge OA 
(Fig. 3.2) and let Oy also lie in this plane. Oz is vertically upwards. At time t, 
let <t> be the angle made by OA with Ox and let 6 be the inclination of the 
plate to Oz. The plate has a component 9 of angular velocity about OA 
and a component tj> about Oz; relative to axes OA, OB, OC, therefore, the 
components of the angular velocity of the plate are (0,<j> cos 9, -<j> s'infl) 
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respectively. Relative to the same axes OABC, the moments and products 
of inertia are 

A = B = jma 2 , C = f ma 2 , 

F = G = 0, H = %ma 2 , 

m being the mass of the plate and a the length of an edge. Using equation 
(3.10), the energy equation for the motion can now be written down: 

6 2 + tf> 2 ( 1 + sin 2 0) - j6<j> cos 9 = 3g( 1 - cos 0)/a. (3.12) 

A second equation is needed to determine the motion. Since the weight 
and the reactions of the pivot and plane have no moment about Oz, 
equation (2.35) requires that h 2 should be conserved. We first calculate the 
components of h in the directions of the axes OABC by using equation (3.1) 
and obtain 

h A = ma 2 (^9—l<f> cos 6), h B = ma 2 (—$9+%<j> cos 0), 
h c = — \ma 2 <j> sin 9. 

Resolving these along the axes Oxyz, leads to the result 
h z = ma 2 [~i<) cos 0+^(1 + sin 2 0)]. 

The equation of conservation of h z is accordingly 

6 cos 0 = j0( 1 + sin 2 0). (3.13) 

Equations (3.12), (3.13), determine the motion. Eliminating <j> between 
these equations, it is found that 

ad 2 (7 + 25 sin 2 0) = 48gr( 1 — cos 0X1 +sin 2 0). 

In principle, this can be integrated to give 9 as a function of t; then, sub¬ 
stituting into equation (3.13), the dependence of 0 upon f follows. □ 

3.3 Motion of Tops 

Problem 3.4 One end of the axis of a top is freely pivoted at a fixed 
point 0 and the top is set spinning. The moment of inertia of the top about 
its axis is C and the moment of inertia about any perpendicular axis 
through 0 is A. Calculate the equation of angular momentum. 

Solution. This problem can be analysed by constructing Euler’s equations 
and the energy equation. We prefer, however, to use it to illustrate the 
application of direct vector methods, to the solution of dynamical problems. 

In Fig. 3.3, Oxyz is a fixed rectangular frame with Oz vertically upwards. 
The unit vector along the axis of the top is a = OA. If to is the vector 
angular velocity of the top about 0, the velocity of A will be to x a. But 
this velocity is also a and hence 

to x a = a. 
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Taking the vector product of both sides of this equation with a and 
expanding the triple product, we get 

a 2 c«> —a.cua = a xa. 

Since a 2 = 1, this yields 

(0 = axa + na, 

where n = a.o> is the component of the top’s angular velocity in the 
direction of the axis. 

Resolving along principal axes at O (i.e. the top’s axis and any two 
perpendicular axes), the angular momentum of the top about O is given by 
h = (A(o i , Ao} 2 , Cn) 

= A(o} l ,(o 2 ,n)+{C-A)(0,0,n) 

= Ao) + (C — A)na 
= Aa x a + CVia. 

If g is the vector acceleration due to gravity, the weight of the top is Mg- 
this acts at the centre of mass C on the axis of the top distant / from 0 
The only other external force acting upon the top is the reaction at the 
pivot 0. Hence, taking moments of these forces about 0 and employing 
equation (2.35), we obtain 

, d 

lax Mg = — (Aaxa + Cna) = Aa x a + Cmi + Cnk. 

(N.B. a x a = 0.) This is the equation of angular momentum. 

Differentiation of the equation a 2 = 1 yields a . a = 0. Thus, if the 
scalar product of both members of the equation of angular momentum is 
taken with a, we find that Cn = 0. (N.B. Both the scalar triple products 
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[aag], [aaa], vanish.) Thus n is constant and the angular momentum 
equation reduces to 

/laxa + Cna = M/axg. (3.14) 

This equation possesses a steady state solution in which the top axis 
rotates with constant angular velocity Si about Oz, making a constant 
angle 9 with the vertical. This motion is termed steady precession. For such 
a motion, the point A describes a horizontal circle of radius sin 9 with 
angular velocity Si. Its velocity a is accordingly Si sin 0 b, where b is a unit 
vector normal to the vertical plane containing the axis of the top (b = OB 
in Fig. 3.3). The acceleration a of A is Si 2 sin 9 along the inwards radius of 
/4’s path, a x a can now be calculated to be the vector — fi 2 sin 9 cos 0b and 
a x g to be g sin 9 b. Thus, equation (3.14) is satisfied provided 9 and Si are 
related by 

-ASi 2 sin9 cos 9 + CnSi sin 9 = Mlg sin 9, 
i.e. ASi 2 cos9-CnSi + Mlg = 0. (3.15) 

Provided C 2 n 2 ^ AAMgl cos 0, this gives two possible angular velocities 
of steady precession for a given inclination 9 of the axis to the vertical and 
a given rate of spin n of the top. □ 

Problem 3.5 Show that, in the general motion of a top, the axis moves 
between two cones. 

Solution. Taking the vector product of both sides of equation (3.14) with 
a, and expanding the triple products, we get * 

/4a. aa = M/a.ga. 

Thus Ai .'i-Mlk . g = 0, 

and this equation integrates immediately to give 

$Az 2 — Min. g = E, 

where E is constant. This is the energy equation. If <f> is the angle made by 
the vertical plane n containing the top’s axis with the plane Oxz, A will 
have velocity components (i) 9 perpendicular to OA and in the plane n, 
(ii) (j> sin 9 in the direction of b. Hence 

a 2 = 6 2 +<j> 2 sm 2 9 

and the energy equation can be written 

$A(d 2 + (j> 2 sin 2 9) + Mgl cos 9 = E. (3.16) 

If k is the unit vector along Oz, since k and g are parallel, taking scalar 
products of both members of equation (3.14) with k, we get 

/4k. a x a + Cnk. a = 0, 

35 














or 


^(/Ik.axa + Cnk.a) = 0. 

Now kx a = sin 0 b and, hence, [kaa] = <£sin 2 0. Thus, integrating the 
last equation, we get 

/t^sin 2 0+Cncos0 = H, (3.17) 

where H is constant. Eliminating (j> between equations (3.16), (3.17), we 
obtain 

At) 1 = 2£-2M g )cosg-L f '- C "'” s9 > i ,, 18) 

Asin 2 8 ' ' ' * 

Putting X — cos 0, this can be written 

X 2 = (a - bX){ 1 - A 2 ) - (c - dX) 2 = f(X), 
where a = 2E/A, b = 2Mgl/A, c = H/A and d = Cn/A. 

f(X) is a cubic in X which must be positive for some value of X in the 
interval (-1,1) (otherwise X is imaginary for all values of 9) Let this 
value be X 0 . Then 

< ftto) > 0> /(l) < 0, /( + oo) > 0. 

P follows that the equation f(X) = 0 has three real roots in the intervals 
( 1» )> (*o > 1)- (1. °o)• If , X 2 are the roots in the first two intervals X 

must always satisfy the inequalities A, < ^ A 2 . Since /. can only vanish 
or X Xj , X - X 2 , the sign of X can only change for these values of X. 
Thus, X will oscillate between these extreme values and 8 will oscillate 
between the corresponding values 0,, d 2 . The behaviour of the top’s axis 
is now clear; it will precess about the vertical Oz with angular velocity 6 
determined by equation (3.17) and, at the same time, it will oscillate 
between two cones with common vertex O, common axis Oz and semi- 
vertical angles 9 X , 0 2 . ^ 

3.4 Rolling and Slipping 

Problem 3.6 A uniform sphere moves on a rough horizontal plane 
Determine its motion (i) when it rolls, (ii) when it slips. 

Solution. Oxyz is a frame fixed in the horizontal plane as shown in Fig. 3 4 
Let ( P, Q, R) be the components of the reaction of the plane on the sphere 
at its point of contact C and let (u, v, 0) be the components of the velocity 
of its centre G. Then, if m is the mass of the sphere, the equations of motion 
of G are 

P = mu, Q = mv, R-mg = 0. (3.19) 

Let (cUj, oj 2 , cuj) be the components of the angular velocity of the 
sphere about G. Since the axes parallel to Oxyz through G are always 
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principal axes for the sphere, its angular momentum about G is given by 
h G = sfno 2 ((o l ,co 2 ,co 3 ), a being the sphere’s radius. Taking moments of 
the external forces about G and employing equation (2.38), it follows that 

fma 2 (ri>, ,(b 2 ,cb 3 ) = (Qa, - Pa, 0). 

Thus, ai 3 is constant, and 

P = — |wiacu 2 , Q = | macu,, (3.20) 

If the sphere is rolling upon the plane, C must be stationary. But. the 
velocity of C relative to G has components ( — aco 2 , aaj,, 0). Since G has 
velocity (u,i>,0), the velocity of C is (u — aa> 2 ,v+aa> l ,0) and the rolling 
conditions are 

u-aco 2 = v+ao) l = 0. (3.21) 

Equations (3.19H3.21) determine the rolling motion. Solving, it will be 
found that ti = v = 0 and P = Q = 0, i.e. the sphere rolls in a straight line 
with constant speed. 

If the sphere is slipping, let 8 be the angle made by the velocity of the 
point C on the sphere with the x-axis. Then 

(p + acu, )/(u - aa> 2 ) = tan 8. (3.22) 

We shall assume that the horizontal (or frictional) component of the 
reaction at C is in the opposite direction to this velocity and is of magnitude 
HR, where g is the coefficient of friction. Hence 

P = -gR cos 8, Q=-gR sin 8. (3.23) 

Equations (3.19), (3.20), (3.22), (3.23), determine the motion. Thus 
v + a<b t lQ/2m v+aw l 

u-aco 2 IP/2m u — ao) 2 

and this implies that 

logfu + aai,) = log(u — aw 2 ) + constant. 
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Hence v+acoi 

u-aco 2 ~ tan0 = con stant, (3.24) 

i.e. the direction of slip never varies. Equations (3.19), (3.23), now give 

“ = ~ MS cos 6, v = -fig sin 6, (3 251 

2 n Sm0, (1} 2 =~cosd. ( 3 . 2 6) 

Thus to,, a) 2 increase or decrease uniformly. Ultimately the rolling 

a " d "" ^ •o r„,7S 

3.S R M „i „8 Axes |„ this SKlta „ we sha|| so|w 

on he motion of a rigid body by calculating the components ofltaar fZ 

axT o r bS'Tr ' r °' a,i " 8 frame f * l >°* «« « always pr?n c ,“a 1 
axes for the body. If x is a vector whose components with respect to F are 

known, ,t ,s necessary to distinguish between Sx/tt the ra^hanlTr 
he vector relative to F and dx/i/t, the rate of change of the tecLfrelSv. 

3 , n ° n rot ^ tln g frame instantaneously coincident with F If is the 
angular veloctty of F, i, is proved in Veaor AlgeFra fbfarder, s e "o n 6 . 4 ' 

dx dx 

dt~ft + QxX ■ (3-27) 

u T it 8 n ^ a " d 

linear and angular moujum rts^xively f o r J’ h L b ^.°’ 
^n 6, T 3J A , wheel of mass m and rad ius a rolls on a rough horizontal 

eel (Fig. 3.5). These axes are principal axes for the wheel• the assn™tf>H 
moments of inertia will be taken to be (mA-,m^, 2 mkT 3 d 

ThenV HIT* 3 ° and radius c ’ with angular velocitv Q 

lhen G has an acceleration cQ 2 towards O and if (R in arA t uJ ‘ , 

and frictional components o, the reaction l 

F = mcQ 2 , R = mg . p 2 8) 
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The frame G123 rotates about a vertical axis with angular velocity Q, 
and hence the components of its angular velocity in this frame are given 
by S2 = (Q cos 9,0,—Q sin 0) (0 is the constant inclination of the wheel 
to the vertical). Let v be the angular velocity of the wheel about its axis for 
motion relative to the frame GI23. Then the net angular velocity of the 
wheel has components o = (Q cos 0, 0,v-Q sin 0). The angular momentum 
of the wheel about G is now given by equations (3.2) to be 

h 0 = mk 2 (Q cos 0,0, 2n), 
where n = v — £2 sin 0. Hence 
dh r , d 

~dt ~ dt ^ C ° S ^ + ^ coss ’ n xmk 2 (Q cos0, 0, 2 n) 

= ntk 2 (0, - (2 2 sin 0 cos 0 - 2nQ cos 0, 0). 

Taking moments of the external forces about G2, equation (2.38) leads to 
the result i 

Ra sin 0 — Fa cos 0 = mk 2 Q cos 0(Q sin 0 + 2n). (3.29) 

In the frame G123, C is the point with position vector ( — a, 0,0). Its 
velocity relative to G is, therefore, ox a = (0, - an, 0). G has velocity 
(0, cQ, 0) and the net velocity of C is accordingly (0, cQ - an, 0). But C is 
stationary if the wheel does not slip and hence 

cQ = an (3.30) 

is the rolling condition. Eliminating F, R and n from equation (3.29) by 
the use of equations (3.28), (3.30), the following relationship between Q 
and 0 is derived, 

aa 2 tan0 

(3.31) □ 


Q 2 = 


k 2 (a sin 0 + 2c) + a 2 c' 
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Problem 3.8 A gyroscope is mounted in gimbals so that its axis is free to 
rotate in a horizontal plane. Show that it behaves like a compass. 

Solution. GIEN are axes through the mass centre of the flywheel such 
that G1 is vertically upwards and GE, GN lie in a horizontal plane pointing 
due east and north respectively. Let G3 be the axis of the flywheel and let 
G2 be taken horizontally to complete the frame G123 (Fig. 3.6). Let A be 



the latitude of G and 0 the inclination of G3 to GN. The frame G123 has 
angular velocities (i) ft about an axis parallel to the earth s axis due to the 
earth’s rotation and (ii) -0 about G1 due to the motion of the gyroscope’s 
axis. Resolving these angular velocities along the axes of the frame, we 
obtained the components (ft sin A — 0, — ft cos A sin 0, ft cos A cos 6) for the 
net angular velocity of the frame G123. 

Let ij/ be the angle through which the flywheel has rotated with respect 
to the frame G123. Then, the angular velocity of the flywheel has com¬ 
ponents (ft sinA—0, - ft cos A sin 0, + ft cos A cos 9) and its angular mo¬ 
mentum about G is given by 

h c = [A(ft sin A — 0), - Aft cos A sin 0, Cn], 
where A, A, C are the principal moments of inertia in the frame G123 and 
n = \j/ + ft cos A cos 9. 

Relative to a non-rotating frame instantaneously coincident with G123, 
the rate of change of h c is 

^ = — [A(ftsinA-0), — Aft cos A sin 0, Cn] 
dt dt L 

+ (ftsinA-0, - ft cos A sin 0, ft cos A cos 0) 
x [ A{Q sin A — 0), - Aft cos A sin 0, Cn] 

= (-Ad, - Aft cos A cos 00, Cn) 

+ [ - Cnft cos A sin 0 + Aft 2 cos 2 A sin 0 cos 0, 

(Aft cos A cos 0 — CnXft sin A — 6), 0]. 
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Thus the equations of angular momentum are 

A0 + Cnft cos A sin 0-Aft 2 cos 2 A sin 0 cos 0 = 0, 

- 2Aft cos A cos 00 +Aft 2 cos A sin A cos 0-Cn(ft sin A-0) = M, 

Cn = 0, 

where M is the moment about G2 of the forces exerted by the gimbals 
(these forces have no moments about G1 or G3). 

The first of these equations can be approximated to read 

A0+Cnft cos A sin 0 = 0, (3.32) 

since ft 2 will be small by comparison with nft. The third equation shows 
that n is constant. 

Equation (3.32) indicates that the axis of the gyroscope will oscillate in 
a horizontal plane about the north line 0 = 0 and that the period of small 
oscillations is 2n^J(A/CnQ cos A). These oscillations will be damped by 
friction and the axis will ultimately come to rest pointing due north. In 
practice, the frictional couple about G3 tending to reduce the spin n is 
cancelled by the torque of an electric motor. □ 

Problem 3.9 A sphere of mass m, radius b and radius of gyration k about 
any diameter, rolls, without slipping, on the outside of a fixed sphere of 
radius a. Determine the possible motions. 

3 



Solution. Take axes G123 through the centre of the moving sphere, 
where G3 is along the common radius OG of the two spheres, G1 lies in 
the vertical plane containing this radius and G2 is horizontal (Fig. 3.7). 
9 is the angle between OG and the upward vertical and 4> is the angle made 
by the vertical plane G13 with a fixed vertical plane. (F, ,F 2 ,R) are the 
components of the reaction of the fixed sphere on the moving sphere at 
the point of contact C. to = (to, ,to 2 ,to 3 ) is the angular velocity of the 
moving sphere. 
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The velocity of G is given by 

v = (a+bW, <i> sin 0,0). (3.33) 

Since C is the point b = (0,0, -b) in the frame G123, its velocity is 
v+ruxb = [(a + i>)$ — bco 2 , (a + fc)0sin0+iKO,, 0]. 

If there is no slipping, we accordingly have the rolling conditions 

bw l = —(a + b)4> sin0, ba> 2 = (a + b)d. (3.34) 

The frame G123 has simultaneous angular velocities (i) 6 about G2 and 
(ii) (j> about a vertical axis. Thus resolving along the axes of the frame 

(2 = (—<t> sm0,6,(j> cos 9). (3.35) 

The acceleration of G in an inertial frame can now be calculated, thus: 

d\ dv _ 

— = — + ft xv 
dt dt 

= {a + b)(9 , $ sin 9 + 6<i> cos 9, 0) 

+(a + b){-(j) 2 sia9cos9, 6<j> cos 9, -6 2 -<f> 2 sin 2 9). 

Resolving the external forces upon the sphere in the directions of the axes, 
the equations of motion of G are now found to be 

m(a + b)(9 — <j> 2 sin 9 cos 9) = F l +mgsin9, (3.36) 

m(a+b)($ sin 9+26<j> cos 9) = F 2 , (3.37) 

— m(a+b)(d 2 + (j) 2 sin 2 9) = R — mgcos9. (3.38) 


Since the axes of the frame G123 are principal axes for the sphere, its 
angular momentum about G is given by 

h G = mfc 2 (co,, co 2 , co 3 ). (3.39) 

Taking moments of the external forces about G, the equation of angular 
momentum is 

(dh G /dt) + ftxh G = (F 2 b, -F,b,0). 

This leads to the equations 

mk 2 (cb,+a) 3 9—(o 2 <fi cos 9) = F 2 b, (3.40) 

mk 2 ((b 2 +o} l <j>cos9+a) 3 (j>sin9) = — Fjb, (3.41) 

mk 2 ((b 3 — cj 2 <$> sin 9— to, 9) = 0. (3.42) 

Equations (3.34), (3.36) to (3.38), (3.40) to (3.42), determine the motion. 
Equations (3.34), (3.42), show that u) 3 = n (a constant). 

Eliminating F,, F 2 between equations (3.36), (3.37), (3.40), (3.41), and 
employing the rolling conditions (3.34), we get 

(k 2 + b 2 )(a+ h)($ sin 9+29<p cos 9) = k 2 hn6, (3.43) 
(k 2 + b 2 )(a + b)(d - (j> 2 sin 9 cos 0) + k 2 bn<j> sin 9 = gb 2 sin 0. (3.44) 
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Multiplying equation (3.43) by sin 9, it integrates immediately to yield 

(k 2 + b 2 )(a + b)<j>sin 2 9+k 2 bncos9 = H, (3.45) 

where H is constant. Multiplying equation (3.43) by <^sin0 and equation 
(3.44) by 9 and adding, we obtain an equation which also integrates 
immediately to give 

%k 2 +b 2 )(a+b)(d 2 + <j> 2 sin 2 9)+gb 2 cos 9 = E, (3.46) 

where £ is the constant of integration. The reader may verify that this is 
the energy equation. 

If equations (3.45), (3.46), are compared with equations (3.17), (3.16), 
respectively, for the motion of a top, it will be observed that, with 
A = m(k 2 + b 2 ), C = mk 2 b/(a + b), I = b 2 /(a + b), they become identical. 
The common radius OG therefore always moves like the axis of a top, 
i.e. it moves between two cones with common vertex O and common 
axis the vertical through 0. This means that the point of contact C moves 
between a pair of horizontal small circles on the fixed sphere. 

In the particular case of steady motion, 9 and <j> are constant and equa¬ 
tions (3.43), (3.44) are satisfied provided 

(k 2 + b 2 )(a+b)<i> 2 cos 9-k 2 bruj>+b 2 g = 0. 

This equation possesses real roots in <j> if 

n 2 > 4g(k 2 + b 2 )(a + b)cos9/k 4 . □ 

EXERCISES 

1. A rigid body whose principal moments of inertia are A, A,C(C > A) 
rotates about its centre of mass under the action of a retarding couple 
- K(o, where K is constant and <o is its angular velocity. The initial values 
of the angular velocity components a>,, a> 2 , <u 3 (relative to the principal 
axes) are (f2 l ,0,f2 3 ) respectively. Show that, at a time t later, 

cu, = Q, e - *' M cosA, co 2 = e”*' M sinA, m 3 = f2 3 e"*' /c , 
where A = — ^~^ Q 3 (1 -e~ K " c ). 

2. A point O is a perpendicular distance c from a smooth plane vertical 
wall. A uniform rod of length a has one end smoothly pivoted at O and 
the other end moves on the plane wall. If 9 is the inclination of the rod to 
the horizontal and <f> is the angle between the vertical plane containing the 
rod and the vertical plane perpendicular to the wall, show that 











a cos 6 cos (f) = c and that the energy equation can be written 


. f csin 2 0 1 3 g . 

6 2 \ 1 4— 77—2 -IZ— T\ = —(sm a- sin 0), 

J{a 2 cos 2 0-c 2 )J a 


y/(a 2 cos 

where initially 0 = a, 6 = 0. 


3. By differentiating equation (3.18), obtain the following equation for 
the motion of a top: 

A('d- (j> 2 sin 0 cos 9)+ Cn<j> sin 0 = M glsmO. 

Putting 0 = a + c, <j> = Si + q, in this equation and equation (3.17) and 
neglecting second order terms in s, q, obtain equations governing small 
oscillations about the steady state determined by equation (3.15) (with 
0 = a). Hence show that the frequency of these oscillations is p/2n, where 

A 2 p 2 = A 2 Q 2 - 2AMgl cos ot + (M 2 g 2 l 2 /Q 2 ). 


4. If in Problem 3.6, the horizontal plane is constrained to rotate about 
Oz with angular velocity Si, but the frame Oxyz remains fixed, show that 
the equations of linear and angular momentum for the sphere remain 
unchanged, but (assuming no slipping) the rolling conditions become 

u-aa) 2 = — Siy, v + ao) l = Six, 

where (x,y, z) are the coordinates of the centre G of the sphere. Deduce that 
G describes an ellipse, making one circuit in a time 7n/Si. 


5. A uniform solid sphere of radius a rolls without slipping inside a fixed 
circular cylinder of radius b, whose axis is horizontal. If 0 is the inclination 
to the vertical of the plane containing the axis of the cylinder and the centre 
G of the sphere, u is the component of G’s velocity along the cylinder’s 
axis and a> 3 is the component of the sphere’s angular velocity along CG, 
where C is the point of contact of the sphere and cylinder, obtained the 
equations 


7(6 — a)6 2 = lOg cos 0 +const., 


lu 2 + 2a 2 o)\ = const., 


d 2 cj 3 2 

- 0 


(Hint: Choose axes G123, where G1 is parallel to the cylinder’s axis and 
G3 is along CG.) 
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Chapter 4 

Equations of Lagrange and Hamilton 


4.1 Lagrange’s Equations If the configuration of a dynamical system 
can be specified at any given time t by stating the values of n independent 
quantities q x ,q 2 ,...,q„, these are said to be generalized coordinates for 
the system. Thus, for the system described in Problem 2.7 (Fig. 2.6), the 
distance x moved by the ring at time t from its initial position and the 
angle 0 made by the rod with the vertical, are suitable generalized co¬ 
ordinates. 

The kinetic energy of the system T can always be expressed as a function 
of the q t and their time derivatives q t (the generalized velocities). The 
dependence of T upon the q t can take any form, but the kinetic energy is 
always a quadratic function of the q { . T may also depend explicitly upon 
t if the system’s constraints (see below) are variable. Thus, in general, 

T = T (<h .9.»4i»—» 4»t) (4.1) 

The particles of a system are not usually free to execute purely arbitrary 
motions. These motions are required to satisfy certain geometrical con¬ 
ditions called constraints. For example, the two ends of the rod in Fig. 2.3 
are constrained to move, one along a vertical line, and the other along a 
horizontal line; further, since the rod is rigid, any two of its particles must 
always remain at the same distance from one another. In cases where 
these geometrical conditions do not change with t, the constraints are 
said to be fixed ; in cases where they do depend upon t, the constraints are 
said to be variable, e.g. if in Fig. 2.3 the vertical wall were to be moved in a 
predetermined manner, this constraint would be variable. If all the con¬ 
straints are fixed, T will be a homogeneous quadratic function of the q t . 

Lagrange’s equations of motion for a system whose kinetic energy has 
been found in the form (4.1) are 



i = \,2,...,n, where Q, are the generalized components of force. The Q, 
are calculated as follows. Suppose the system to be in an arbitrary configur¬ 
ation (q x ,q 2 ,...,q n ) at time t and let it undergo a small displacement 
to a neighbouring configuration (q, +Sq l ,q 2 + dq 2 ,...,q n +Sq n ). This 
displacement must be consistent with the constraints, assumed frozen in 
in the form they take at the instant t. Thus, if the constraints are variable, 
this displacement cannot coincide with the actual motion during the interval 
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(t, f+<5t); for this reason, the displacement is referred to as a virtual displace¬ 
ment. During the virtual displacement, all forces of the system (internal and 
external) are supposed to remain constant in magnitude and direction. 
The work done by these forces is always expressible, to the first order, in 
the form 


W = t Q,Sq t . (4.3) 

i= 1 

The coefficient of Sq, in this expression accordingly defines the component 
of force Q t . 

It should be noted that the work done by the internal forces of a rigid 
body during any virtual displacement which does not disrupt the body, is 
zero. Similarly, provided the length of a string is not altered in a virtual 
displacement, its internal tension does no work. 

Problem 4.1 Let (r, 0,0) be the usual spherical polar coordinates of a 
particle m with respect to a rectangular inertial frame. Let (F r ,F 9 ,F + ) be 
the components of the force acting upon the particle in the directions r 
increasing, 6 increasing and 0 increasing, respectively. Find Lagrange’s 
equations for the motion. 

Solution. In a virtual displacement, the particle’s coordinates change to 
(r+<5r,0+<50,0+<50) and the work done by the force is given by 

<51 V = F r <5r+F # r 50+ F^r sin 050. 

It follows that 


Qr = F r , Qo = rF e , Q, = rsin0F*. 

The spherical polar components of the velocity of the particle are (r, 
r0, r sin 0 0) and the particle’s kinetic energy is accordingly given by 

T = T(r,0,r) = jm(r 2 + r 2 6 2 + r 2 sin 2 0 <j> 2 ). 


We can now write down Lagrange’s equations associated with the 
generalized coordinates (r, 0,0), 

m(r -rd 2 -r sin 2 00 2 ) = F r , 
m^-^(r 2 0)—rsin0cos0<^ 2 J = F e , 


m 


rsin0 dt 


(r 2 sin 2 00) = F. 


The expressions multiplying m in the left-hand members of these equations 
are, of course, the spherical polar components of the acceleration vector. 

□ 
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Problem 4.2 A light string passes over a pulley of mass 4 m which can 
rotate freely about a fixed horizontal axis. A particle of mass 2m is attached 
to one end and a pulley of mass m to the other. The second pulley is also 
free to rotate about a horizontal axis and carries a light string with masses 
2m, 11m at its two ends. The system moves in a vertical plane with all 
strings vertical. Calculate the accelerations of the masses and pulleys. 




Solution. At tinje t, suppose the particle A (Fig. 4.1) has moved upwards a 
distance x and the particles B, C have each moved a distance y relative 
to the lower pulley. Then x, y are generalized coordinates for the system. 

Consider a virtual displacement in which x increases by Sx and y by by. 
The work done by the weight of A is then - 2mg Sx and by the weight of the 
lower pulley is mg Sx. B moves down a distance Sx—Sy and its weight does 
work 2 mg(Sx—Sy). C moves down 5x + 5y and its weight does work 
Umg(Sx + Sy). None of the remaining forces of the system do work. Thus 

SW = I2mg Sx + 9mg Sy, 
whence Q x = 12 mg, Q y = 9 mg. 

The angular velocity of the lower pulley is y/b, where b is its radius; 
the velocity of its centre of mass is x. Assuming it to be a uniform disc, 
its kinetic energy is accordingly |mx 2 +±my 2 . The kinetic energy of the other 
pulley is mx 2 . The particles B, C have velocities (x — y), (x + y), respectively 
in the downwards direction. Thus the kinetic energy of the system is 

T = £mx 2 +£my 2 +mx 2 + m(x-y) 2 +Y»>(;c+y) 2 + mx 2 
= |m(4x 2 + 4xy + 3y 2 ) 

= T(x,y). 
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Lagrange's equations can now be written down in the form, 

^(18x + 9j>) = 12 g, J t (9x + ~y) = 9 g. 

Whence, x = \g, y = \g, from which the accelerations of all parts of the 
system can be calculated. q 

Problem 4.3 For the system described in Problem 2.7, take x, the distance 
moved by the ring and 9, the rod s inclination to the vertical, as generalized 
coordinates. 

Solution. 

T =T{9, x, 6) = jmx 2 + m(x 2 + a 2 6 2 - 2a6x cos 9) + $ma 2 d 2 . 

In the virtual displacement, x increases by Sx and 9 by 69. Only the 
weight of the rod does work and, since its centre of mass moves a distance 
a 69 sin 9 vertically, hence 

6W = - 2mga 69 sin 9. 

Thus Q x = 0, Q 0 = — 2mga sin 9. 

Lagrange’s equations can now be constructed: 

— (3x—2 ad cos 9) = 0, — a 2 6 —2ax cos 9^j - 2 adx sin 9 = — 2ga sin 9. 

Since x = v, the first of these equations is equivalent to the equation of 
conservation of horizontal momentum (2.26). The second equation can 
be written 

I a6-x cos9+gsin9 = 0. 

After putting x = \ad cos 9, this is found to be identical with the result of 
differentiating equation (2.28). g 

4.2 Systems with Variable Constraints 

Problem 4.4 A rigid straight rail is fixed in the plane of a lamina and a 
small engine, of mass m, moves itself by friction drive along the rail. The 



Fig. 4.2 
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lamina is horizontal and is free to rotate about a vertical axis which does 
not intersect the rail. Discuss the motion of the system. 

Solution. In Fig. 4.2, £ is the engine and EN is the rail; 0 is the point in 
which the vertical axis of rotation intersects the lamina. ON = a is 
perpendicular to the rail and EN = x. We shall assume that the motion 
of the engine relative to the rail is known, i.e. x is a known function of t; 
this constitutes an internal constraint of the system. 

Let <t> be the angle turned through by the lamina; this is the generalized 
coordinate for the system. If mk 2 is the moment of inertia of the lamina 
and rail about the axis through 0, its kinetic energy is %mk 2 <f> 2 . The engine’s 
velocity is compounded of its velocity x along NE relative to the lamina 
and the velocity of the lamina at £, namely r(f> perpendicular to OE 
(r = OE). Thus, the net kinetic energy of the system is given by 

T = jmk 2 <f> 2 + %m(x 2 + r 2 <j) 2 — 2r<f)xcos9), 

where 9 = angle NOE. Since r cos 9 — a and r 2 = x 2 +a 2 , this reduces to 

T = %m(x 2 + a 2 + k 2 )<j> 2 - max<j>+%mx 2 . 


In a virtual displacement, 0 increases by 8<l> and the engine remains 
stationary on its rail. Thus the engine’s motor does no work and neither 
do any ofthe other forces belonging to the system. Hence = 0. Lagrange’s 
equation for the system is now found to be 

^[(x 2 + a 2 + k 2 )(f> - ax] = 0. 

Hence (x 2 + a 2 + k 2 )<j> = ax. 


This equation determines (f>(t) when x(r) is given. It may, however, be 
integrated, without any knowledge of x(t),thus: 


4 > 


-l 


a dx 


tan 


-1 


x 2 + a 2 + k 2~ + * 2 ) J {a 2 + - 


where we have assumed that <f> = 0 when x = 0. This equation relates 
the motions of the engine and lamina. In particular, as x varies from — oo 
to +oo, it is clear that <j> varies from —na/2j(a 2 + k 2 )lo +na/2 s /(a 2 + k 2 ) 
and, hence, that the lamina can never rotate through an angle exceeding 
na/J(a 2 + k 2 ). □ 


4.3 Conservative Systems If the system is conservative, a potential 
energy function V exists. V will be expressible in terms of the coordinates 
< 7 , and the time t may occur explicitly if the constraints are variable. The 
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work done by the forces of the system in a virtual displacement in which 
g t increases by Sq, is given by 


5 W= -SV = - 


Hence 

and Lagrange’s equations are 


0 „J_1 

Q ‘ V 


i — 


dV 

dq, 


5q t . 


d_ fST\_dT _ _3V 
di \dq,J dq, dq,' 


Defining a Lagrangian function L by the equation 

L = L(q i ,'..,q„,q l ,...,q n ,t) = T-V, 
Lagrange’s equations (4.2) can be expressed in the form 


d_ (^\_8L _ 

dt \dq,J dq, ~ 


(4.4) 

(4.5) 


(4.6) 


Problem 4.5 Derive the energy equation from Lagrange’s equations 
tn the case of a conservative system with fixed constraints. 

Solution. For a conservative system with fixed constraints, Tis a homo¬ 
geneous quadratic function of the generalized velocities. Thus 

27 = £ (4.7) 

where the a,j are functions of the coordinates. It follows that 

V . dT 

L*'d$= 2T (4.8) 

Multiplying equation (4.5) by q, and summing with respect to i. 


or 

since V = 



V ( Q i.42»- -»4«)- Now 


(4.9) 


T = 7X9it-..9 a .4i. 4 „\ 

if the constraints are not variable (i.e. T is not explicitly dependent on t) and 
thus 


dT Y f . dT _ dTl 
dt 
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Using this equation and equation (4.8), equation (4.9) can accordingly be 
written 


2 — — — - - — 
dt dt ~ dt 

and this is equivalent to the energy equation T+ V = constant. □ 


Problem 4.6 A uniform rod is smoothly pivoted at one end to a fixed 
point and moves three dimensionally under gravity. Obtain its equations 
of motion. 



Solution. Let 0 be the inclination of the rod to the vertical and 4> the 
angle made by the vertical plane containing the rod with a fixed vertical 
plane. Take axes OA along the rod, OB horizontally and OC to complete a 
rectangular frame. These are principal axes for the rod at the pivot 0 
and its principal moments of inertia are 0, fma 2 , f ma 2 , 2 a being the 
length of the rod. Its angular velocity about 0 has components (4> cos 0, 
6,(j>smQ) and its kinetic energy is accordingly given by 

2 T = %ma 2 (d 2 + <f> 2 s\n 2 6). 


V = — mga cos 0 and the Lagrangian is thus 

L = %ma 2 (d 2 + <j> 2 sm 2 9)+mgacos0. 
Lagrange’s equations are 

’()-<£ 2 sin0cos0+^sin0 = 0, — (<j>sin 2 6) = 0. 
4a dt 


Integrating the second equation, we get 

4> sin 2 0 = h, 

which also follows from the fact that angular momentum about the vertical 
through 0 is conserved. The first Lagrange equation can be replaced by the 
simpler energy equation; however, for the steady motion in which 6 = a 
(a constant), this equation yields the important result 


t 2 


3g 

— sec a, 
4a 


determining the angular rate <j> in terms of the inclination a. 


□ 
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Problem 4.7 Find the equations of motion for the two particle system 
studied in Problem 2.1. 

Solution, r and 0 are generalized coordinates. The potential energy of the 
hanging particle can be taken to be - kmg(l — r), where I is the length of 
the string. Thus, neglecting a constant term, 

V = kmgr. 

Also T = jm(r 2 + r 2 6 2 ) + %kmr 2 . 

Thus L = £m(fc +1 )r 2 + \mr 2 9 2 - kmgr. 

Lagrange’s equations are 

[m(/c +1 )r] - mr9 2 + kmg = 0, ~(mr 2 9) = 0. 

The first of these equations follows from equations (2.1), (2.3) by addition 
and the second is equation (2.2). □ 

Problem 4.8 FindtheequationsofmotionforthetopstudiedinProblems 
3.4 and 3.5. 

Solution. Generalized coordinates for the top are 0, <p, where 0 and <p 
are shown in Fig. 3.3 and ^ is the angle between a plane fixed in the top 
containing the top’s axis OA and the vertical plane OAN (also indicated 
in the figure). The rectangular frame OABC has angular velocity compo¬ 
nents (<j> cos 0,0, <j> sin 0) in the directions of its axes and the top has 
angular velocity (i£,0,0) relative to this frame. The angular velocity of the 
top relative to a fixed frame instantaneously coincident with OABC is 
accordingly {if> + <f> cos 0,0, 4> sin 0). Since the axes of this frame are principal 
axes of inertia for the top, its kinetic energy is now calculated to be 
$A(6 2 + (j> 2 sin 2 0) + $C(<j> cos 0+ 1 l/) 2 . 

The potential energy of the top is Mgl cos 0. Hence 

L = %A(0 2 + <£ 2 sin 2 0) + \C(<i> cos 0 + \j/) 2 - Mgl cos 0 
and Lagrange’s equations take the form 

— (Ad) —A<i> 2 sin 9 cos 9 +C(<j> cos 9 +sin 9—Mgl sin 9 = 0. 
dt 

^-[A<j>sin 2 9 + C(<j>cos9 + \j/)cos9'\ = 0, ^-[C($cos0+i£)] = 0. 

dt dt 

The last of these equations integrates to give 

<j>cos9 + ii/ = n, (4.10) 

where, as in Problem 3.4, n is the constant component of the top’s angular 

52 


velocity in the direction of its axis. The second Lagrange equation then 
integrates to yield 

A<f>sin 2 9+Cncos9 = H, (4.11) 

This is equation (3.17). The first Lagrange equation can be written 

A’O—A<j) 2 sin 9 cos 9+Cn<j> sin 9 = Mgl sin 0. (4.12) 

Usually, it is convenient to replace this equation by the energy equation, 
which has the advantage that it is of first order instead of second order. 
However, for steady precessional motion 0 = a of the top, equation (4.12) 
reduces to the important equation (3.15) immediately. □ 

Problem 4.9 Find the equations of motion of the falling square plate of 
Problem 3.3. 

Solution. The kinetic energy of the plate has already been calculated. 
Its potential energy is %mga cos 0. Hence 

L = %ma 2 [d 2 + 4> 2 {l + sin 2 0) -10<£> cos 0] - \mga cos 0. 

0, (f> are, of course, generalized coordinates for this system. Lagrange’s 
equations are now calculated in the form 

— (d-i<j>cos9) — (j> 2 sin9cos9— f0$sin0 = —sin0, 
dt a 

+sin 2 0)-|^cos0] = 0. 
dt 

The second of these equations is equivalent to equation (3.13). The first 
equation is best replaced by the simpler energy equation, (3.12). □ 

If the Lagrangian does not contain one of the generalized coordinates 
explicitly, e.g. q t , then the associated Lagrange equation simplifies to the 
form 



and is immediately integrable to yield 

dLjdq l = constant. (4.13) 

dL/dq x is called the component of generalized momentum of the system 
associated with the coordinate q x and equation (4.13) can therefore be 
interpreted as a law of conservation of momentum. As an example, the 
Lagrangian for a top is independent of both <t> and 4 1 This has the conse¬ 
quence that the associated Lagrange equations integrate to yield equations 
(4.10), (4.11). It is left for the reader to verify that these equations imply 
that the components of the angular momentum of the top in the directions 
of the top’s axis and of the vertical are conserved. 
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4.4 Relaxation of Constraints If a particular force belonging to a system 
does no work during a virtual displacement, it will not enter into any of 
the Lagrange equations and hence cannot be calculated by this method. 
Such a force is the reaction P of the smooth wall in problem 2.3. If it is 
desired to calculate such a force employing the Lagrange technique, 
it is necessary to imagine that the associated constraint is abolished and 
that the force is supplied by an alternative external agent. As a consequence, 
the number of coordinates is increased. Equations of motion for the new 
system are then written down and these will now include the required 
force. The manner in which this force must then be chosen in order that 
the system’s motion satisfies the original constraint can now be found. 
This technique is termed the method of relaxation of constraints. 

Problem 4.10 Calculate the reaction P of the wall on the rod for the 
system described'in Problem 2.3 (Fig. 2.3). 

Solution. Remove the wall. Relaxation of this constraint means that the 
horizontal motion of the rod is not determined by its rotational motion 
and the point 0 beneath the end A of the rod can move horizontally. 
Let x be the distance of 0 to the left of a fixed point on the horizontal 
ground at any time t. Then (x, 6) are generalized coordinates for the relaxed 
system. Let x be increased by Sx and 6 by S6 in a virtual displacement. 
The work done is 6 W = MgaSOsind+PSx.ThusQ e = Mga sin 0,Q X = P. 

The centre of mass of the rod has.velocity a() perpendicular to OG 
relative to 0. This must be added vectorially to the horizontal velocity 
x of O. The kinetic energy of the rod is therefore given by 

T = iMaW +1 M(x 2 + a 2 d 2 + 2adx cos 9). 

Lagrange’s equations are accordingly 

fa'd+x cos 9 = g sin 0, M(x + ad cos 9 - ad 2 sin 9) = P. 

We now choose P so that x is constant. The above equations then become 
ja'd = g sin 0, P = M(ad cos 9— a9 2 sin 9). 

Since d = d(jd 2 )/d9, the first equation can be integrated to yield the energy 
equation 

f ad 2 = g(cosa-cos9) 

and P now follows exactly as in Problem 2.5. □ 

Problem 4.11 Calculate the force exerted on the particle C by the 
string supporting it in the system of Fig. 4.1. 

Solution. We imagine the string cut at C and equal, but opposite, forces 
1 applied to C and to the cut end of the string. Let y denote the distance B 
has moved upwards relative to the lower pulley and let z denote the 
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distance C has moved downwards with respect to this pulley. Having 
relaxed the constraint, y and z are now independent coordinates. The 
kinetic energy is found to be given by 

T = m[9x 2 +%y 2 +*fz 2 -2xy+ llxz]. 

In a virtual displacement for which x, y, z increase by Sx, by, bz respec¬ 
tively, the work done by the two forces F is F(Sx + <5y) ~F(bx + bz) = 
F{by-bz). Thus, the net work for the whole system is 

bW = F(by-bz)-2mgbx + mgbx + 2mg(bx-Sy)+Umg(bx+bz) 
and it follows that 

Q x = 12 mg, Q y = F-2mg, Q z = -F+llmg. 

Lagrange’s equations are 

18x—2j)+ llz = \2g, h-2x = --2g, Uz + 1 lx = -- +1 lg. 

2 m m 

Putting z = y, these equations are easily solved to give 
X = 20 > y = 19, F = *fmg. 

The first two of these equations have already been obtained in Problem 
4.2. The third equation giving the force applied to C bythe string, can be 
obtained very simply from the fact that C’s acceleration is known to be 
.x+y = %g; its equation of motion is accordingly 

\lmg-F =\lmfy). □ 

4.5 Non-Holonomic Systems The mechanical systems studied in the 
previous sections of this chapter have all possessed the property that the 
increments bq,(i = 1 , 2 ,..., n) in the coordinates determining a virtual 
displacement could be chosen arbitrarily. Such systems are said to be 
holonomic. This property is crucial for the argument by which Lagrange’s 
equations are established and these equations are accordingly valid for 
holonomic systems only. 

Consider a circular disc of radius a constrained to roll without slipping 
on a horizontal plane with its plane vertical. Taking axes Oxy parallel to 
the horizontal plane, let (x,y) be the coordinates of the centre of the disc. 
Let 9 be the angle made by the plane of the disc with the x-axis and let <p 
be the angle through which the disc has rotated. Then (x, y, 9, <p) constitute 
a set of generalized coordinates for the disc. If, however, <f> is increased by 
b<f>, the centre of the disc moves a distance a S(f> in a direction making an 
angle 9 with Ox and the increments in x and y are determined by the 
equations 

Sx = a Sep cos 9, by = a 6<p sin 0. (4.14) 
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This means that the increments Sx , 5y, 8<f> are not independent of one 
another and, hence, cannot be chosen to take arbitrary values. Neverthe¬ 
less, it is not possible to integrate the equations (4.14) to yield relations 
between the four coordinates (x, y, 0, <t>), which could be used to reduce the 
number of coordinates necessary to specify a configuration of the disc on 
the plane. For this system, therefore, although there are four coordinates, 
the increments of only two of these (<50 and 6<f>) can be chosen arbitrarily 
in a virtual displacement. 

In the general case, if q t , q 2 ,..., q n are the generalized coordinates and 
the virtual increments Sq t , 8q 2 ,...,8q„ are related by r (< n) independent 
conditions, then n — r can be chosen arbitrarily and the system is said to 
have n-r degrees of freedom. For a holonomic system, r = 0 and the 
number of degrees of freedom is equal to the number of generalized 
coordinates. 

Let the r conditions satisfied by the 6q ( be 


X ^tj — 0> 


(4.15) 


(= i 


where j = 1,2 ,...,r and the A tJ are functions of the q,. Then, Lagrange’s 
equations for the system are 

df8T\ dT „ ' A , ...» 

■ e,+ £ ( 1 

where i = 1,2,_, m. The r quantities Xj are called Lagrange's multipliers ; 

they are functions of t which are eliminated from the calculation by making 
use of the r relations between the generalized velocities q ( which follow 
from equations (4.15). For example, if the constraints are fixed, the virtual 
displacement is a possible actual displacement of the system and hence, 
dividing the equations (4.15) by St , we obtain 


X A,jq t = 0 


(4.17) 


i= 1 


as a set of conditions on the generalized velocities. 


Problem 4.12 A wheel, of mass m and radius a, rolls without slipping on 
a horizontal plane. Obtain its equations of motion. 

Solution. Oxy are fixed axes in the horizontal plane (Fig. 4.4). The plane 
of the wheel intersects the horizontal plane in a line / making an angle 0 
with Ox. 0 is the inclination of the wheel to the vertical. Cl 23 are rectangular 
axes through the centre of the wheel such that Cl passes through its point 
of contact P with the horizontal plane, C2 is parallel to / and C3 is per¬ 
pendicular to the plane of the wheel. Relative to these axes, the wheel has 
rotated through an angle 0. Then, if (x,y,z) are the coordinates of C 
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relative to the axes Oxyz , the five quantities (x,y, 0,0,0) are generalized 
coordinates for the wheel. 

The frame Cl23 has an angular velocity d about C2 and an angular 
velocity 0 about the vertical through C. Thus, the net angular velocity of 
the frame has components (<j> cos 0, d y 0 sin0) along its axes. The wheel 
has angular velocity 0 about C3 relative to this frame. Hence, the net 
angular velocity of the wheel has components (0 cos 0,0, 0 + 0 sin 0). 
The position vector of P relative to C has components ( — a, 0,0). By 
taking the vector product of these two vectors, we find the velocity of P 
relative to C, (0, - a(0 + 0 sin 0), ad). Resolving these components in the 
directions of the axes Oxyz, we get new components 

- ad cos 0 sin 0 - a(0 + 0 sin 0)cos 0 along Ox, 

a6 cos 0 cos 0 - a(0 + (j> sin 0)sin 0 along Oy, 

ad sin 0 along Oz. 

But C has velocity (x,y,z) in the frame Oxyz. The components of the net 
velocity of P in the frame Oxyz can now be found and these must all 
vanish. Thus, writing i^ + 0sin0 = n, 

x —ndcos0sin 0 —nncos0 = 0, (4.18) 

y + a0cos0cos0 — an sin 0 = 0, (4.19) 

z+n0sin0 = 0. (4.20) 

These are the rolling conditions at P . Since z = acosO , equation (4.20) is 
obvious. 

57 

E 


















Since the constraints are fixed, a virtual displacement can be an actual 
displacement and the first two rolling conditions show that the virtual 
increments 8x , etc., are related by the conditions 


Sx — aS9 cos 0 sin 0 — a(<50 + 8(f> sin 9 )cos 0 = 0, 


+ a<50 cos 0 cos 0 — a(<50 + <50 sin 0)sin 0 = 0. 

These are the forms taken by equations (4.15) for this problem; the matrix 
with elements A {j is therefore as set out below: 


1 0 
0 1 


Mu) 


- a cos 9 sin 0 
a sin 9 cos 0 
a cos 0 


a cos 0 cos 0 
a sin 0 sin 0 
a sin 0 _ 


(4.21) 


C123 are principal axes for the wheel. Taking the corresponding 
principal moments of inertia to be (mk 2 ,mk 2 ,2m/c 2 ), the kinetic energy of 
the rotary motion of the wheel about C can be found using formula (3.1 1 ). 
The Lagrangian then follows in the form 

L = fm/c 2 [0 2 cos 2 0 + 0 2 + 2(0 + 0 sin 0) 2 ] 

+ %m(x 2 + y 2 + a 2 0 2 sin 2 0) — mga cos 0. (4.22) 

(Note: Use has been made of equation 4.20.) 

Lagrange’s equations (4.16) associated with the coordinates (x, y, 0,0,0) 
are respectively 

^(mx) = A f \ j t (my) = k 2 , (4.23) 

■7 [m(fc 2 +a 2 sin 2 0 )f)] — mk 2 <f> cos 0 ( 2 n — 0 sin 0 ) 
dt 

— ma 2 0 2 sin0 cos 6 —mga sin 0 = — a cos 0(A, sin 0—A 2 cos 0), (4.24) 

^[mfc 2 (<£cos 2 0 + 2nsin0)] = — asin0(A, cos 0 + A . 2 sin 0), (4.25) 

2 mk 2 n = — a(A, cos 0 + A 2 sin 0). (4.26) 

The first two of these equations show that the multipliers A,, A 2 can be 
interpreted, physically, as the x- and y-components of the plane’s reaction 
upon the wheel at P. The rather complex third equation may, conveniently, 
be replaced by the energy equation. These equations, together with the 
rolling conditions (4.18), (4.19), determine the motion of the wheel. 

The general case of the motion of the wheel will not be studied further 


58 


here (see, however, Exercise 5 at the end of this chapter). Instead, we shall 
consider two special cases, (i) small oscillations about steady rolling along 
a straight line and (ii) steady rolling around a circle. 

For steady rolling along a straight line, the quantities x, y, 0 , ri, 0, A,, x 2 , 
all vanish. It is easy to verify that this constitutes a special solution of the 
equations of motion. If small oscillations are taking place about this 
motion, all these quantities can be regarded as small. To the first order of 
small quantities, Lagrange’s equations approximate to the forms 

mx = A,, my = A 2 , 

mk 2 9- 2mk 2 n<j> - mgaO = - a(A , sin 0 - X 2 cos 0), 

0 + 2nd = 0 , 2 mk 2 n — —a(X i cos 0 + k 2 sin 0 ). 

Differentiating the rolling conditions (4.18), (4.19), and retaining first 
order terms only, we find that 

x = aO sin 0 + an cos 0 — an<p sin 0 , 
y = —aid cos (f> + an sin 0 + an<j> cos 0 . 

Eliminating , A 2 , x, y, the following equations are obtained: 

(k 2 + a 2 )d—(2k 2 + a 2 )n<j> = ga6, 0 + 2 n 0 = 0 , n = 0 . 

Thus n , the angular velocity of the wheel about its axis, can be treated as a 
constant. Integrating the second equation to give 0 + 2n0 = 0 (assume 
0 = 0 when 0 = 0 ) and substituting for 0 in the first equation, we get 

d+p 2 0 = 0 , 

where p> . 

V k 2 +a 2 

The condition for small oscillations is now seen to be 

n 2 > ga/2(2k 2 + a 2 ). 

In the case of a hoop, k 2 = \a 2 and the condition is n 2 > g/4a. For a 
uniform disc, k 2 = $a 2 and the condition is n 2 > g/3a. 

For steady motion in a circle 0, <f> and n are all constant. Thus, the 
rolling conditions reduce to 

x = an cos 0 , y = an sin 0 . 

Differentiating, x = — an 0 sin 0 , y = an<j> cos 0 . 

A,, A 2 can now be calculated and, then, the equation (4.24) gives 

k 2 0 cos 0 ( 2 n — <^sin 0 )+ 0 asin 0 = — a 2 ntj> cos 0. (4.27) 
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Suppose the centre of the wheel describes a circle of radius c with angular 
velocity Q. Then (j> = -Q (in Fig. 4.5, $ will be decreasing) and an = cQ. 
Equation (4.27) accordingly reduces to 

ga 2 tan 0 

<c 2 ( 2 c+qsin 0)+a 2 c 

which has already been derived as equation (3.31). □ 

4.6 Hamilton’s Equations For any conservative system the n quantities 
dL/dq, are called the generalized components of momentum and are denoted 
by p,. Thus 

p, = dL/dq,. (4.28) 

Since all terms in L are of the second degree at most in the q t , p, depends 
linearly upon these generalized velocities. It is therefore possible to solve 
for the q, as linear expressions in the p, , thus: 

4i = '£ B'jpj+Ci, (4.29) 

J m i 

B tj , C, being functions of the coordinates q, and t. 

Evidently the configuration of the system and the velocities of all its 
parts, i.e. the dynamical state of the system, can be specified by stating the 
values of the q, and p ,. q, and p, are said to be conjugate quantities. 

The Hamiltonian of the system is defined by the equation 

H = £ Pj4(-T, (4.30) 

i= j 

where it is assumed that all the velocities q, are replaced by linear expres¬ 
sions in the momenta using equation (4.29). Thus 

H = H(p, , q ( , t). (4.31) 

Problem 4.13 Show that the following equations are equivalent to 
Lagrange’s equations 

4 , = dH/dp„ p, = -dHJdq,. (4.32) 

Solution. From equation (4.30), treating the 4/ as functions of the basic 
independent variables (p,, q, , t) (equation 4.29), it follows that 
dH _ y dq, dL y dL dq, 
dq } ~ Y'dq, dqj 
Using equation (4.28), this reduces to 


dH 

dL 

d (dL 

dq } ~ 

1 

CDl 

C? 1 
II 

dt\dqj 
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where the last equality follows from Lagrange's equations. The second of 
equations (4.32) has now been proved. 

Again from equation (4.30), we have 


dH . v V dL dq, 

to-*'* L^-L^-1r 


up, ->• y>,- 

This is the first of equations (4.32). Equations (4.32) are called Hamilton's 
canonical equations of motion. □ 


Problem 4.14 Show that, if a system has fixed constraints, then H is the 
total energy. 

Solution. In such a case, the kinetic energy is a homogeneous quadratic 
form in the q t (equation (4.7)). It follows that 

P, = dL/dq, = dT/dq, = £ a tj qj. 

J 

Substituting in equation (4.30), this gives 

H = Y i a tJ q,4 J -L = 2T-(T-V) = T+V. n 

Problem 4.15 Find Hamilton’s equations for a top. 

Solution. The Lagrangian for a top has been given in Problem 4.8. The 
generalized components of momentum are derived from this in the forms 
p 9 = dL/dd = Ad, 

Pq = dL/d(j> = A(j) sin 2 0 + C((p cos 0 + t/<)cos 0, 
p^ = dL/d\j> = C(<j> cos 0 +if/). 

Solving for the generalized velocities, we obtain the equations 


9 = pJA 

(4.33) 

<i> = (P 4 -P 4 , cos 6 )/(A s\n 2 0) 

(4.34) 

, cos 6 [ 1 cos 2 0 1 


^ ~ 4sin 2 fl P * + L c + ^ sin2e J P *’ 

(4.35) 


Since the constraints are fixed, the Hamiltonian can be identified with 
the total energy. Thus 

H = $A(d 2 + <j> 2 s\n 2 0)+$C(<j>cos0 + il/) 2 + Mgl cos 0 

= £4 Cp« + (p^ - p* cos O) 2 cosec 2 0 ] + ^pl + Mgl cos 6. 

Hamilton’s equations can now be constructed. The first three are simply 
the equations (4.33H4.35); the remainder are as follows: 

Pe = ^ (P* ~ P<, cos cos 0-p^ )coscc 3 9 + Mgl sin 0, 


P* = P* = 0 . 
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The first of these equations is easily seen to be equivalent to the Lagrange 
equation (4.12); the remaining equations integrate to yield equations which 
are equivalent to equations (4.11) and (4.10). □ 

Problem 4.16 A bead of mass m is threaded on a smooth circular wire 
of radius a and the wire is constrained to rotate in its own plane, which is 
horizontal, about a point 0 of itself with constant angular velocity m. 
Discuss the motion of the bead. 

Solution. OA is the diameter of the wire through the pivot 0 and P is the 
position of the bead at time t. If C is the centre of the wire, angle 
PC A = 2 x angle POC = 26. The velocity of the bead relative to a frame 
rotating with the wire is 2a& perpendicular to CP. The velocity of this 
frame at P is to.OP = 2am cos 6 perpendicular to OP. The velocity of 
the bead in a fixed frame is found by adding these velocities vectorially; 
thus the kinetic energy of the bead is 

T = 2ma 2 (6 2 +m 2 cos 2 0 + 2mdcos8). 

Since the bead’s potential energy is constant, T is also the Lagrangian 
for the bead. As 6 is the only coordinate, the conjugate momentum is 

p = 8T/d6 = 4mu 2 (()+«cos 0). 



The constraint is variable and H must therefore be calculated from 
equation (4.30), 

H = p6—2ma 2 (6 2 +m 2 cos 2 8 + 2m6cos0) = (p 2 /ima 2 )-pm cos 6. 
Hamilton's equations are 

& = (p/4ma 2 )—cocos 0, p = —pm sin 6. 

Eliminating p, we calculate that 

20 + co 2 sin20 = 0, (4.36) 
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proving that the angle PCA oscillates about zero in synchrony with a 
simple pendulum of length g/m 2 . n 


Problem 4.17 If F = F{p l ,q t ,t) is a quantity whose value depends on 
the system’s state, prove that 


dF 

dt 




where 
Solution. We have 


[f ’ H]= 1 (t-t-itt)- 

M \ dc h dPi fy, dqj 


dF _ 8F_ . dF . dF 
dt ~ dj l qi + di, p,+ ~di 

and the result now follows immediately by use of Hamilton’s equations. 

□ 

[F,H] is called a Poisson bracket. In particular, if F = H, the Poisson 
bracket vanishes and we have dH/dt = 8H/3t. Thus, if H is not explicitly 
dependent upon t, it is a constant of the motion. In particular, if the con¬ 
straints are fixed, this implies that the energy of the system is conserved. 
However, in the previous example, the energy of the bead is not conserved 
and yet 

H = g~ 2 “ P“>cos9 = 2ma 2 (6 2 —m 2 cos 2 6) = constant. 

The last equation also follows by integration of equation (4.36). 


EXERCISES 

1. A uniform rod of mass m and length 2 a has a small light ring attached 
to one end. The ring is threaded on to a smooth vertical rigid wire and the 
rod moves in a vertical plane under gravity. If, at time t, the ring is a 
distance x below a fixed point on the wire and 6 is the angle made in the 
downwards direction by the rod with the horizontal through the ring, 
obtain Lagrange’s equations and deduce that 

x = gt -ad cos 6+ constant, (4-3 cos 2 0)6 2 = constant. 

2. Four uniform rods, each of mass m and length 2a, are smoothly 
joined together at their ends to form a rhombus ABCD. AB is smoothly 
pivoted at its centre and the configuration moves in a vertical plane with 
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DC below AB. If 9 is the angle made by AB and DC with the horizontal 
and <f> is the angle made by AD and BC with the vertical, show that the 
Lagrangian for the system is given by 

L = faui 2 (2& 2 + 50 2 )+ 4mga cos <(>. 

Deduce that AB and DC rotate with constant angular velocity and that 
AD and BC oscillate in synchrony with a simple pendulum of length 5a/3. 

3. The ends of a uniform rod of mass 3m and length 2a are constrained to 
move along two smooth wires OA and OB, OA being horizontal and OB 
vertically downwards. A smooth bead of mass m is free to slide on the rod. 
The whole system is made to rotate about OB with constant angular 
velocity yj(6g/5a). If 9 is the inclination of the rod to OA, obtain Lagrange’s 
equations and show that these possess a solution corresponding to steady 
motion with the bead at rest at the mid-point of the wire and sin 0 = f 


a fixed vertical diameter. A smooth bead of mass m is threaded on the wire. 
Show that the Hamiltonian for the system is 

H = 

where <j> is the angle through which the wire has rotated and 9 is the angle 
made by the radius to the bead with the downward vertical. Obtain 
Hamilton s equations and derive from these the equations 

0(1 +sin 2 0) = constant, 0-0 2 sin0cos0 = -gsinO/a. 

Show that a steady motion is possible with <j> = a, 9 = a and acu 2 = gseca. 


4. AC, BC are two equal uniform bars, smoothly hinged at C. A light 
ring attached to A is threaded on a smooth horizontal wire Ox and a 
similar ring attached to B moves on a smooth vertical wire Oy, pointing 
upwards. The bars move in the vertical plane Oxy with C above Ox. If 
9, (f> are the respective inclinations of AC, BC to the horizontal (measured 
in the same sense), calculate Lagrange’s equations satisfied by these 
coordinates. If, initially, 9 = fa, (f> = fa, show that, at this instant, the 
angular accelerations of the bars are 9^/ 11a and 3 y /‘3g/22a. 

5. In the problem of the wheel rolling on a horizontal plane (Problem 
4.12), eliminate x, y, 2,, A 2 from the Lagrange equations (4.23H4.26) and 
the rolling conditions (4.1*8), (4.19), and hence obtain the equations 

(k 2 + a 2 )Q - (2k 2 + a 2 )n(j> cos 9 + k 2 <j> 2 sin 9 cos 9 = gasinfl, 

0 cos 9 — 200 sin 9 + 2nd = 0, (2k 2 + a 2 )n + a 2 6<i> cos 9 = 0. 


Deduce that n satisfies the equation 

where n = sin 9. 

6. Approximate the equations for a rolling wheel obtained in the 
previous exercise in the case where 0 and n are small and 0 = co + c, 
where co is constant and e is small. Hence show that a wheel spinning with 
angular velocity cu about a vertical diameter is in stable motion provided 
co 2 > ga/(k 2 +a 2 ). 

7. A uniform circular wire of radius a and mass 2m is free to rotate about 
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Chapter 5 

Small Oscillations about Equilibrium 

5.1 Stability of Equilibrium Consider a holonomic conservative system 
whose constraints are fixed. If there exists a configuration in which it can 
rest in equilibrium, then its equations of motion (4.5) must possess a 
solution in which the coordinates < 7 , are all constant and the velocities < 7 , 
all vanish. Since the constraints are fixed, T will be given by equation (4.7) 
and it follows that 

dT/dq, = dT/eq, = 0 

identically for this particular solution. Substituting in equation (4.5), we 
deduce that 

dV/dqf = 0, (5-D 

i = l,2,...,n, in the position of equilibrium; i.e. V is stationary in such a 
position. 

If the position of equilibrium is stable, then V is a relative minimum in 
this position; for otherwise it will be possible for the system, after a small 
disturbance, to move away from the position in such a way as to decrease 
V and, by the energy equation, this will lead to an increase in T and thus 
to a more rapid movement of the system away from equilibrium. 

Problem 5.1 A uniform rod AB of length 2 a and mass m has its end A 
connected to a point C of a smooth vertical wall by an elastic string of 
natural length \a^jl and Hooke’s constant mg/a. The end B of the rod 
rests against the wall vertically below C. Show that there is only one 
equilibrium configuration and that it is unstable. 
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Solution. Consider a configuration in which B is a distance d below C 
and the length of the string is /. Then, if 0 is the angle indicated in Fig. 5 . 1 , 
by the cosine rule 

cost) = (l 2 -4a 2 -d 2 )/4ad. ( 5 . 2 ) 

The potential energy of a string whose Hooke’s constant is k and whose 
extension is x, is jkx~. Thus the potential energy of our configuration is 
given by 

V = - mg(d+a cos 0)+^ (l-\oy/1) 2 


= -^(3^ + / 2 -4a 2 )+g(/-W7) 2 , 

having substituted for cos 0 from equation ( 5 . 2 ). 

We now introduce dimensionless variables x, y by putting d = ax and 
/ = ay. Then 

V = —3x + 2(y —| v / 7 ) 2 J. ( 5 . 3 ) 

For equilibrium, V must be stationary with respect to variations of x 
and y. The conditions for this are 

V-4 ,1 n 

—2 - 3 J = °’ (5.4) 

•~ + 4 0'-iV 7 )J = 0 . (5.5) 

Equation (5.5) shows that 

y = Xyjlf(2x - \ ). ( 5 . 6 ) 

Substituting from this equation for y into equation ( 5 . 4 ), we obtain 
3x 4 —3x 3 + 3x 2 —4x+l = 0 . 

By inspection, x = 1 is seen to be a root and taking out the factor (x - 1) 

we get h 

(x-lX3x 3 + 3x-l) = 0. 

If ./M = 3x +3x—1, then /'(x) = 9x 2 + 3 > Oand hence /"(x) increases 
monotomcally. It follows that /(x) can have but one zero and, since 

= s* tflis lies in the interval (0,£). But for values of x in 
this lnterva! equation (5.6) shows that y is negative and this is impossible 
Thus the only configuration for which V is stationary is the one deter¬ 
mined by x = Uy = yj~l, '\.e.d = a, l = a^7,0 = ^7t. In this configuration 
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= < mga 

dV . 
Ty =Vnga 












Thus 


d 2 V/8x 2 = mga(4-y 2 )/2x 2 = - 3mga/2, 
d 2 V/6y 2 = mga(2x- l)/2x = mqa/2, 
d 2 V/dxdy = mgay/2x 2 = Jlmga/2. 


d 2 VT (PVcPV 
dxdy\ dx 2 dy 2 


= f mga > 0 


and the stationary value of V is neither a maximum nor a minimum. It 
follows that the equilibrium is unstable. □ 


Problem 5.2 Equal point electric charges are fixed at the vertices of a 
square. A similar charge is constrained to move in the plane of the square. 
Show that it can rest in stable equilibrium at the centre of the square. 

Solution. Taking axes Oxy through the square’s centre parallel to its 
sides, the coordinates of the vertices can be taken to be (±a, ±a). Let 
(x,y) be the coordinates of the movable charge when it is displaced 
slightly from 0. The potential energy of a movable charge e when at a 
distance r from an equal fixed charge is e 2 /iie 2 /4ne 0 if S.I. units are used). 
Thus the potential energy of e at (x, y) in the field of the equal charge at 
(a, a) is 


e 2 [(x-a ) 2 + (y-a) 2 ] * = 

e 2 

J2a 


±-A*‘+y‘)] * 


l--(x + y)+ 
a 

l+^(x + y)+g^(x 2 + y 2 + 6xy) 


>} 


to the second order in the small quantities x and y. Three similar expres¬ 
sions can be derived for the contributions of the charges at the other three 
vertices. The net potential energy of the movable charge is calculated to 
be given by 


V = 


2 v /2 a 


3 ( 8 fl 2 + x 2 + y 2 ). 


It is now evident that V is a minimum at 0 and hence that the movable 
charge is in stable equilibrium at this point. 


5.2 Small Oscillations about Stable Equilibrium If a holonomic con¬ 
servative system with fixed constraints is slightly disturbed from a position 
of stable equilibrium, it will commence to oscillate about this position. 
We shall choose the generalized coordinates q , in such a way that these 
are all zero in the equilibrium position and shall assume that, together 
with the velocities q,, they remain small during the oscillatory motion. 
By neglecting terms of order higher than the second in the Lagrangian, an 
approximate solution for the oscillatory motion can be found. 
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Problem 53 A light string is stretched between two fixed points A, D 
on a smooth horizontal table with tension P. Two particles B , C of mass m 
are fixed to the string at the points of trisection. Discuss possible oscillatory 
motions of the particles in the plane of the table, at right angles to the string. 



Solution. Let x, y be displacements of the particles from their equilibrium 
positions (Fig. 5.2). If x, y are small, the tensions in the three segments of 
the string will remain P to the first order. Consider the component in the 
direction of the displacement x of the tension in BC acting upon the 
particle B; this is P sin a = Ptanot = P(y-x)/a, where/4Z) = 3a. Similarly, 
the same component of the tension in BA is -Px/a. The net force acting 
upon B in the direction of the displacement x is accordingly P(y—2x)/a. 
Similarly, the net force acting upon C in the direction of y is P(x-2y)/a. 
Hence, if x increases by Sx and y by <5y, the work done by the forces of the 
system is 

P[(y-2x)6x+(x-2y)dy]/a 

and it follows that the generalized forces for the system are given by 
Qx = P(y-2x)/a, Q y = P(x-2y)/a. 

The kinetic energy of the system is given by 

T = im(x 2 +y 2 ) 

and Lagrange’s equations are therefore 

-£ ( mx) = £(y—2x), l(my) = £(x-2y). 

Putting a > 2 = P/ma, these can be written 

x + w 2 (2x—y) = y+tu 2 (2y—x) = 0. (5.7) 

It is physically evident that these equations possess an oscillatory solution 
and we accordingly seek a solution in the form 

x = /4sin(/»+a), y = Bsin(pt+a). 

Upon substitution, it is found that the equations are satisfied provided 
(p 2 — 2co 2 )A + to 2 B = 0, ( 5 . 8 ) 

co 2 A+(p 2 —2 co 2 )B = 0. ( 5 . 9 ) 
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For these equations to have solutions other than A = B = 0 (the equili¬ 
brium solution), it is necessary that 

p 2 — 2w 2 to 2 _ q 

(O 2 p 2 — 2 (D 2 

i.e. (p 2 - co 2 Xp 2 — 3eo 2 ) = 0. (5.10) 

Thus, p = wot .J3 (o (negative values of p do not provide distinct solutions). 
In the case p = to, either of the equations (5.8), (5.9), show that A = B 

and the solution found is 

x = y = A sin(cot+a), (5.11) 

A and a being arbitrary. The configuration of the particles in this type of 
motion is shown in Fig. 5.3(a); the particles oscillate in phase with one 
another with period In/cD. 



In the case p = v /3 (o, equations (5.8), (5.9), show that A = - B and the 
motion is determined by 

x = — y = A sinQlwt + ai). (5.12) 

This type of motion is illustrated in Fig. 5.3(b); the particles oscillate in 
anti-phase with one another with period In/Jiw. 

These two possible types of oscillatory motion are called the normal 
modes of vibration of the system. The characteristic feature of a normal 
mode is that all particles of the system oscillate in phase or in antiphase 
with one another. 

Since the equations (5.7) governing the system’s behaviour are linear, it 
follows that the solutions (5.11), (5.12) may be added to yield a more 
general solution, 

x = A sin(cut-t-a)-l-/4'sin( v / 3a>H-a'),) ^ ^ 

y = /lsin(cut + a)-/4'sin( N /3«t-l-a').J 

This solution contains four arbitrary constants A, A', a, a' and is the general 
solution of the equations (5.7). Given the initial positions of the particles 
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and their velocities, these four constants can be found and the subsequent 
motion determined. For example, if the system is started from rest at 
t = 0 with x = a,y = -3a, then A = a, A' = 2a, a = 37t/2,a' = 7t/2,i.e. 

x = a( — cos ojt + 2 cosV 3a>l), y = — a(cos a>f +2 cos ^ 3a»t). □ 
Problem 5.4 Two uniform rods having the same length 2a and masses 
2m, 3 m, are smoothly hinged together at one end of each. The free end of 
the less massive rod is smoothly hinged to a fixed point, so that the rods 
are free to swing in a vertical plane under gravity. Calculate the periods of 
the normal modes of small oscillations. 



Solution. Let 9, <f> be the inclinations of the rods AB, BC to the vertical 
at time t (Fig. 5.4). The velocity of the pivot B is 2 aO and the velocity of 
the centre of mass of BC relative to B is a<j>. The net velocity v of the centre 
of mass of BC is obtained by summing these velocities vectorially we 
find that 

v 2 = a 2 [40 2 + <j> 2 +40<p cos {<t>— 0)]. 

An expression for the kinetic energy of the system can now be constructed, 
T = jma 2 0 2 +$ma 2 <j> 2 +%mv 2 . 

Substituting for v 2 and subtracting the potential energy, we generate 
the Lagrangian, 

L = ma 2 [2ld 2 /3 + 20 2 + 6d<j> cos(4> - 0 )] 

+ 2 mga cos 6+3mga(2 cos 0 + cos <f>). 

We next approximate this expression to the second order of small quan¬ 
tities 9, <j>, 6, <j>. The result is 

L = ma 2 ( + 2<^ 2 + 6d<j>)— mga(49 2 + \<j> 2 ). (5.14) 
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Lagrange’s equations for the oscillations are now found to be 

^0 + 6 $+—0 = 0, 60 + 4$ + —0 = 0. 

3 a a 

Taking a normal mode solution in the form 

0 = A sin(pf+ot), $ = Bsin(pf + a), 

we find after substitution that A, B must satisfy the equations 

~ — 4p 2 Jb = 0. (5.15) 

For a non-null solution, it is necessary that 

Hence p 2 = 3g/a or 63 / 17 a and the periods of the normal modes are 
accordingly 

2n^J{a/3g), 2nJ(l 7a/6g). 

If p 2 = 3g/a = o>\, equations (5.15) have the solution A:B = 1:—2 
and the normal mode oscillation is given by 

0 = asin(a>, f+a), $ = — 2 asin(£o, t + a). 

If p 2 = 6g/\7a = a> 2 , then A:B = 3:4 and the normal mode is 
0 = 3asin(a > 2 t+a), $ = 4a sin(a > 2 1 + a). 

The general solution is obtained by superposition of these normal 
modes, 

0 = asin(a), t + a) + 3a'sin(co 2 r+a'), ) 

<f> = — 2asin(cu, t+a)+4a'sin(tu 2 r+a'). J 

□ 

5.3 Normal Coordinates 

Problem S.S Calculate the normal coordinates of the double pendulum 
studied in the previous problem. 

Solution. Referring to equations (5.16), it is evident that if we introduce a 
new pair of generalized coordinates x, y for the system, related to the old 
pair 0 , (f> by the linear transformation equations 

0 = x + 3 y, <t> = -2x+4 y, (5.17) 

then the general oscillatory motion of the system is expressed by the equa¬ 
tions 

x = asin(<u, f + a), y = a' sin(ta 2 1 + a'). (5.18) 

In particular, in the first normal mode, a' = 0 and, hence, y is identically 
zero; in the second normal mode a = 0 and x is identically zero. 


-6p 2 B = 0, —6p 2 /4 + j^ 
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Further, if we substitute from equations (5.17) into the Lagrangian (5.14), 
it reduces to the form 

L = ^ma 2 (x 2 + 51 y 2 ) -10 mga{x 2 + 6y 2 ), (5.19) 

i.e. both quadratic forms representing the kinetic and potential energies of 
the system have been reduced to a sum of squares ( diagonal form). 

The simplification can be carried a little further by introducing coordi¬ 
nates «, v, related to x, y by the equations 

w = V(¥ ma2 *)> v = v /(340wa 2 .v). (5.20) 

In terms of u and v, the Lagrangian (5.19) then assumes the form 

L = 5 <m 2 + u 2 )— 0 (f « 2 + -frv 2 )/a = f(!i 2 + v 2 - cofu 2 - co|u 2 ). (5.21) 

Thus, if the coordinates u, v are used, the kinetic and potential energies are 
expressed in the canonical forms 

T = i(u 2 + i) 2 ), V = j(co 2 u 2 +ojjv 2 ). (5.22) 

Coordinates which permit the kinetic and potential energies of an 
oscillatory system to be written in these forms are termed normal co¬ 
ordinates. The corresponding Lagrange equations are clearly 

ii + (i) 2 u = 0, v+colv = 0. (5.23) 

One normal mode is now specified by the equations 

u = Psin^it+a), v = 0, (5.24) 

and the other normal mode by the equations 

u = 0, v = P'sin(aj 2 t+a'). (5.25) 

This illustrates another characteristic of normal coordinates, that in a 
normal mode, all but one remain zero. 

To summarize, for the double pendulum, the normal coordinates are 
related to the angles 0, <p by the equations 

0 N /(34Oma 2 ) = i+/(51)+3t;, </> v /(340ma 2 ) = -2w v /(51)+4u, 
i.e. u = (A9—3(f))yJ(ma 2 /l5), v = (20 + $)V(17ma 2 /5). (5.26) □ 

In general, approximating to the second order of small quantities 
g,, q„ the kinetic and potential energies of a system can be written in the 
forms 

^ = *1 Mi4i» V - i £ KWj’ (5-27) 

where a iJt b xj are constants. Both these quadratic forms will be positive 
definite , i.e. will be greater than zero for all sets of values of the q i ,q i (except 
that T = V = 0 for the null set). It is proved in textbooks of algebra that 
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two such positive definite quadratic forms can be reduced to the diagonal 
forms 

T = ±2>f, V = cofuf, (5.28) 

i i 

by a real regular linear transformation (from the coordinates q ( to new 
coordinates w,) of the form 

u i = (5.29) 

j 

where the c i} are constants. u { are the normal coordinates for the system. 
With T and V in the forms (5.28), Lagrange’s equations for the system are 

M,+ft)f«, = 0, (5.30) 

i = 1,2,.... n. Substituting the normal mode solution 

u, = /l,sin(pt+a), (5.31) 

we obtain the equations 

A i (cof-P 2 ) = 0 , (5.32) 

i = l,2,...,n.Thus,p = (o t orAf = O.Ifwetakep = (o i ,ihenA 2 ,A 3 ,...,A n 
must all vanish (we here assume that the <y, are all different). Thus, the first 
normal mode is 

«, = /4, sinfitf! t+a), u 2 = u 3 = ... = u n = 0. (5.33) 

Similarly, in the second normal mode only u 2 is non-zero, and so on. We 
conclude that, in any normal mode, all the normal coordinates vanish 
except one which oscillates sinusoidally. 

Problem 5.6 Two equal uniform straight beams AB, BC, of mass m 
and length 2a, are smoothly hinged together at B. They are supported at 
A, B and C by light vertical springs such that, in equilibrium, A, B, C are 
at the same level. The stiffness of the spring supporting B is double that 
of the springs at A and C. Calculate the normal modes of vertical oscilla¬ 
tion of the beams and determine the normal coordinates of the system. 

Solution. Let k be the Hooke’s constant for the springs at A and C, and 
2k the constant of the spring at B. In equilibrium, the springs at A, C will 
be subjected to compressive forces jmg and will each contract a distance 
mg/2k. The central spring will support a weight mg and will also contract a 
distance mg/2k. At some instant t during motion, let A, B, C have moved 
distances x, y, z vertically upwards respectively from their equilibrium 
positions. The potential energies of the springs in this configuration will 
be jk(x-mg/2k) 2 , k(y-mg/2k) 2 , \k(z-mg/2k) 2 , respectively. The centres 
of the beams will have been raised distances Kx+y), ^(y-l-z) above their 
equilibrium positions. It follows that the potential energies of the weights 
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are \mg(x 4-y), jmg(y+z). The net potential energy of the system is accord¬ 
ingly. 


+> 0 (* + 2 3 ' + z >. 


= 2 k(x 2 + 2y 2 + z 2 )+constant. 


(5.34) 


The velocities of the centres of mass of the beams are j(x+y), 5 (y+i). 
The velocity of A relative to B is x —y and it follows that the angular 
velocity of A B is (x— y)/2a. Similarly, the angular velocity of BC is (y— z)/2a. 
We can now write down the kinetic energy of the beams, 

T = £m(x+y) 2 +£m(y+ z) 2 + &n(x - y) 2 + £ m(y - z) 2 , 

= £m(x 2 + 2y 2 -( -z 2 +xy+yz). (5.35) 


Forming Lagrange’s equations, we get 

2x+y+co 2 x = 0, x + 4y + z+2co 2 y = 0, y+2z+co 2 z = 0, 
where or = 6 k/m. Substituting the normal mode solution 

x = 4sin(pt + oc), y = Bsin(pH-a), z = Csin(pr+a), 
the equations are satisfied provided 


(o) 2 -2p 2 )A- p 2 B = 0, 

— p 2 A +(2w 2 —4p 2 )B—p 2 C = 0, 
-p 2 B + (co 2 -2p 2 )C = 0. J 
For a non-null solution, it is necessary that 

co 2 — 2p 2 —p 2 0 


— p 2 2o) 2 — 4p 2 —p 2 

0 —p 2 (o 2 —2p 2 


= 0 . 


(5.36) 


Thus, p = (o/yj 3, to/yj2 or co. If p = cu/^3, equations (5.36) show that 
A:B:C = 1:1:1 and the normal mode is given by 

x = y = z = asin[(cut/ v /3)+a]. 


In this motion, the beams continue in a horizontal straight line as they 
oscillate. If p = cu/^/2, then A:B:C = 1:0: — 1 and the normal mode is 

x = -z = 6sin[(cof/. s /2)+/S], y = 0. 

In this motion, the points A, C oscillate with equal amplitudes in anti¬ 
phase and B remains stationary. If p = oi,A:B:C = 1 :- 1:1. The normal 
mode is 

x = -y = z = c sinfiut + y). 


A, B and C now all oscillate with equal amplitudes, but B is in anti-phase 
with A and C. 
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The general solution of the equations of motion is given by 

x = a sin cut+a^+fc sin ^ j-a>r+/^+csin(atf+y), 

y = asin^-^wf + a^ -csin(a>f+y). 

cot + ^j-b sin +f^j+c sin (cot + y). 


z = a sin 


These equations suggest a transformation to new coordinates C, t], £, where 

x = S + r,H, y = €-C, z = 

Substituting from these equations into equations (5.34), (5.35), we find that 
V = $k(4i 2 +2r\ 2 +4C 2 ), T = £m(6c 2 +2? + 2£ 2 ). 

Normal coordinates u, v, w can now be constructed by taking 
u = y/(2m)Z % v = w = \ / (! m )C- 

Then V = %foo 2 u 2 +%co 2 v 2 +co 2 w 2 ), T = %(u 2 + v 2 + w 2 ). 

Thus, the normal coordinates are related to the original coordinates 
x, y, z by the equations 

1 “ 7(^0 ‘"-Va-'+v'Sw). □ 


EXERCISES 

1. Eight equal point electric charges e are fixed at the vertices ( + «, 
+ a, ±a) of a cube. Show that the potential energy of another charge e 
when it is at the point (x,y,z) is given approximately by 

V = 3 [ 1 + i^a i * 3y * z2+ 322x2+3x2 y2 ~ x *-y 4 - z 4 ) J’ 

when x, y , z are small. Deduce that the ninth charge cannot rest in stable 
equilibrium at the centre of the cube. If this charge is constrained to move 
along a diagonal of the cube, show that the centre is a position of stable 
equilibrium. 

2. A uniform bar AB of length 18a and mass 6 m is smoothly hinged at A 
so that it can rotate in a vertical plane. A smooth ring of mass 3m is 
threaded on the bar and is connected to a point C at a height 16a vertically 
above A by a light elastic string. The modulus of the string is 14mg and its 
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natural length is 8a. Show that the system has a configuration of equili¬ 
brium in which the angle BAC is 60° and the ring is distant 6 a from A , 
but that the equilibrium is unstable. 

3. A uniform triangular lamina OAB can turn freely in a vertical plane 
about the point 0 which is fixed. The lamina has mass 3m and is isosceles 
with OA = OB = 2a and the angle AOB is 120°. A particle of mass m 
is suspended from A by a string of length a and an equal particle is suspended 
from B by a string also of length a. Small vibrations take place in the vertical 
plane OAB. Show that the frequencies of the normal modes are in the 
ratios ^3:1:2. If 6 is the angle made by AB with the horizontal and </>, ip 
are the angles made by the strings with the downward vertical, all angles 
being in the same sense, show that the normal coordinates are given by 
«1 = (i»w 2 )*W>-^), U 2 = (£ma 2 ) 4 (80 + tf> + ^), « 3 = (jma 2 )*(0 -<f> - ^). 
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Moment, generalized components of, 53, 
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Orthogonal transformation. 7 
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Relaxation of constraints, 54 
Rolling on a plane, 36 
Rotating axes. 38 
Rotation about fixed axis, 16 
Rotation, instantaneous axis of, 32 
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Satellite, period of, 15 
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Spherical polar components of accelera 
tion, 46 

Steady rolling, small oscillations about, 59 

Top,oscillations about steady motion. 44 
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Velocities, generalized. 45 
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